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. ABSTRACT 


This program is conducted to study the effects of impurities and 
defects on the performance and permanence of silicon solar cells both 
theoretically using computer models and experimentally. 

Numerical solutions are obtained from the exact one dimensional 
transmission line circuit model to study the following effects on the 
terrestrial performance of silicon solar cells: (i) interband Auger 
recombination, (ii) surface recombination at the contact interfaces, 

(iii) enhanced metallic impurity solubility, (iv) diffusion profiles, 
and (v) defect-impurity recombination centers. Thermal recombination 
parameters of titanium impurity in silicon are estimated from recent 
experimental data. Based on these parameters, computer model calcula- 
tions show that titanium concentration must be kept below 6x10^ ^ Ti/cm^ 
in order to achieve 16% AMI efficiency in a p+/n/n+ silicon solar cell 
of 250 pm thick and 1.5 ohm-cm resistivity or 5x10^^ donors/cm^. 

Computer studies of 10 ohm-cm high efficiency p+/n/n+ silicon solar 
cells give the following indications. (i) Interband Auger recombination 
is unimportant in both the diffused emitter and back surface field layers 
even if experimental interband Auger recombination rates used in the 
computer model are in error and too small by a factor of 100. (ii) 
Interface recombination at the front and back surfaces cannot be impor- 
tant unless the interface recombination velocities lie between the 
diffusion velocity, 10® cm/s, and the effective recombination velocity, 
10® cm/s, in the diffused emitter or back surface field layer. (iii) 

The difference in cell performance between the Gaussian and the measured 
boron and phosphorus diffusion profiles is negligible. (iv) Enhanced 
solubility of recombination centers in the diffused emitter will decrease 
the short circuit current, while in the back surface field layer, it 
will decrease the open circuit voltage. 15% reduction can be expected 
if the solubility of the recombination center is 10% of the diffusion 
impurities in these layers. (v) A similar reduction can occur if the 
boron related impurity— vacancy recombination center is present in the 
boron diffused emitter or if the phosphorus-vacancy (E-center) recom- 
bination center is present in the phosphorus diffused back surface field 
layer, (vi) Analysis of the experimental data of Sandia high efficiency 
back surface field cells shows that either a 1% gold or 0.1% phosphorus- 
vacancy recombination center in the back surface field layer can account 
for the observed 16.8% AMI efficiency, 620 mV open circuit voltage and 
30 mA/cm^ short circuit current. The maximum attainable AMI efficiency 
of this cell design is 19.26% if both the emitter and back surface field 
layer recombinations- are negligible and base lifetime is 700 ps. 

Analyses of diffusion length, lifetime and transient capacitance 
data of Ti— doped silicon crystals and silicon solar cells indicate that 
there are three recombination levels. The two lower levels, E^+336 and 
Ej,-470 mV, are donor like and Ti related while the upper level, Ej,-182, 
is a complex center containing possible oxygen, vacancy and carbon. 
Thermal capture and emission rates at these Ti levels are estimated from 
these experimental data. 


- 1 - 



1 . 


INTRODUCTION 


The objective of this program is to determine the effects of impurities 
and defects on the performance and permanence of silicon solar cells. It 
includes theoretical (computer model) and experimental studies of the effects of 
impurities on the properties of silicon intentionally doped with specific 
impurity elements, as well as the effects of these impurities on the impurity 
related energy levels, concentration of these energy levels and recombination- 
generation properties of electrons and holes at these energy levels. 

This second technical report contains extension of the theoretical and 

computer analyses of silicon solar cell performance reported in the first 

technical report.^ New factors considered are: diffusion profiles of the 

emitter and back surface field layers, interband Auger recombination, and band- 

bound thermal recombination localized at the impurity and defect-impurity centers 

whose solubilities are enhanced in the diffused emitter and back surface field 

layers. A detailed analysis of the recombination properties of the Ti impurity 

is also made by combining Northrop ’s diffusion length data, Westinghouse's 

concentration and OCD lifetime data and our preliminary transient capacitance 

data. These recombination parameters are then used to compute the theoretical 

solar cell performance which is then compared with the experimental AMI 

performance of Ti-doped n+/p and p+/n silicon cells without anti-reflection 

coating reported by Westinghouse. Projection of the maximum allowable Ti 

concentration for a given AMI efficiency is computed for the highest efficiency 

cell structure, the anti-reflection coated p+/n/n+ diode. To achieve 16% AMI 

12 3 

efficiency, Ti concentration must be kept below about 6x10 Ti/cm . 

Due to the large number of variables and parameters considered in the 
theoretical models and two types of recombination models used in the high 
efficiency cells and Ti-doped cells, the following format will be used 
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in the presentation of the results. The format within each chapter is similar 
to that used in the first technical report.^ The technical discussion are given 
here in two chapters, 2 and 3. Chapter 2 covers the high efficiency cell which 
is used to analyze the effects of diffusion profiles and Auger recombination 
mechanism. Chapter 3 covers the Ti-doped cells and a projection of maximum 
allowed Ti concentration. Se”ction titles, and starting page numbers are given below- 


2. TECHNICAL DISCUSSION - HIGH EFFICIENCY CELLS 

2.0 Objectives 4 

2.1 Approach 5 

2.1.1 Recombination Model 6 

2.1.2 Computer Model 12 

2.1.3 Device and Material Parameters 15 

2.1.4 Solar and Optical Parameters 17 

2.2 Results on Back Surface Field 21 

2.2.1 Effect of Back Surface Field 22 

2.2.2 Effect of Diffusion Profile 26 

2.3 Results on Interband Auger Recombination 27 

2.3.1 Emitter Layer 29 

2.3.2 Back Surface Field Layer 30 

2.4 Results on Surface Recombination 31 

2.4.1 Emitter Surface Recombination 31 

2.4.2 Back Surface Recombination 37 

2.5 Results on Enhanced Impurity Solubility 38 

2.5.1 Emitter Layer 40 

2.5.2 Back Surface Field Layer. 44 

2.6 Results on Defect-Impurity Recombination Centers 46 

2.6.1 Emitter Layer 50 

2.6.2 Back Surface Field Layer 56 

3. TECHNICAL DISCUSSION - Ti-DOPED CELLS 61 
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3.2 Results on p-Base Cells 73 

3.3 Results on n-Base Cells . 79 

3.4 Projected Maximum Performance 80 
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2 . 


TECHNICAL DISCUSSION - HIGH EFFICIENCY CELLS 


2.0 Obj ectives 

The purpose of these numerical calculations is to determine 
the effects of emitter and back-surface-field (BSF) recombination on 
the. efficiency (EFF) , open circuit voltage (VOC) , short circuit current 
(JSC) , and fill factor (FF) , as well as the position dependent internal 
characteristics of silicon solar cells. The internal characteristics 
studied include the steady-state carrier concentration, recombination 
rates, lifetimes, impurity charge states, electric field, electrostatic 
potential and diffusion, drift and recombination currents. The 
position variations of these internal characteristics are particularly 
valuable for a fundamental understanding of the effects of impurity 
on cell performance.^ 

The study covered in this chapter includes the effects of 
interband Auger recombination in the high-concentration diffused 
emitter and .BSF layers, -the-.enhanced- solubility -of -recombination- 
impurities in these layers, the presence of defect-impurity pairs 
or higher order complexes as electron-hole recombination sites, 
as well as the shape of the majority impurity (boron, phosphorus and 
arsenic) diffusion profile on the recombination currents . 
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2.1 Approach 

The transmission line circuit model will be used in the numerical 
solution to give the cell performance characteristics. A one dimensional model 
is used so that the solutions are the maximum achievable since lateral losses 
and two and three dimensional inhomogeneities are not included. The 
transmission line method is particularly advantageous for device physics 
analysis since the physical mechanisms, such as diffusion, recombination, 
optical generation and boundary conditions at the contacts, are directly 
represented by circuit elements (resistance, capacitance, dependent and 
independent current sources) in the model. Thus, variations of these 
physical mechanisms can be readily made by modification of the circuit 
elements. In addition, the transmission line circuit model can be solved 
numerically by the well-known and highly developed and efficient matrix 
techniques which are available in many scientific computer languages. 

The details of the logical and computational sequence of solutions were 
described in the first technical report, (TR-1, section 2.1.3)^ and will 
not be repeated here. Modifications and improvements xfill be described in 
detail here in the following subsections. 
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2.1,1 Recombination Models 


The coupled gold donor- acceptor recombination center is used as the 
recombination model to study the effects of BSF, interband Auger recombination 
and enhanced solubility in subsections 2.2, 2.3 and 2.4. In section 2.5, a 
single level acceptor from the phosphorus-vacancy pair (PV) is used as the 
recombination model in the phosphorus diffused BSF and a single level donor 
complex, involving vacancies and unidentified Impurities, is used as the 
recombination model in the boron diffused emitter. The recombination parameters 
of these centers are given in the tables belovf. 

Table 2.1 Recombination Parameters at the 
Gold Center in Silicon at 297K.^ 

14 3 

Acceptor Level Lifetime at 10 Au/cm 


*^n0 

= 

.30000E-08 

cm^/sec 

Tno = 3.333 

ps 

Snl 

= 

.88500E+03 

1 /sec 



Cpl 

= 

.13743E-06 

cm^/sec 

Tpi = 72.76 

ns 

'®pb 

= 

V50000E+02- 

1- /sec 



Ej-Ei 

= 

.85710E-01 

eV 



Erp 


Eq - 547.23 

meV 



Donor Level 





‘^pO 

= 

.55948E-07 

cut/ sec 

■’^pO =178.7 

ns 

VI 


.17500E+08 

1 /sec 



^n-l 

= 

. 61200E-07 

cm^/sec 

-163.4 

ns 

%i0 

= 

.2I000E-01 

1 /sec 





-.26400E-00 

eV 



Ej' 


Ey + 355.66 

meV 
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Table 2.2 Recombination Parameters at the 
Donor Level in Boron Diffused 
Layer in Silicon at 292K. 


Donor Level 


o 

o 

II 

.18194E-07 

3/ 

cm /sec 

e , = 

.23205E+06 

1 /sec 

p-j. 

c , = 

.11700E-05 

cm /sec 

n^i. 

(D 

O 

II 

.98460E+01 

1 /sec 


-.18211E+00 

eV 


Ey + 354 

meV 


lA 3 

Lifetime at 10 Levels/cm 


T _ = 549.45 ns 

pO 


T , = 8.547 ns 

n-1 


Table 2.3 Recombination Parameters at the 
Phosphorus-Vacancy (PV) Acceptor 
Level in Phosphorus Diffused Layer 
in Silicon at 29 7K. 


Acceptor Level 



Lifetime at 

14 , 

10 PV/cm' 

'^nO 

= 

.27660E-06 

3 / 

cm /sec 

'^nO 

36.153 ns 

®nl 

= 

. 78330E+05 

1 /sec 




= 

.20000E-05 

3 / 

cm /sec 

^Pl 

5.000 ns 

%o 

= 

. 75320E+03 

1 /sec 




= 

.84827E-01 

eV 






E^ - 471.6 

meV 



, T 


C 





The data given in the above two tables are obtained from voltage stimulated 
transient capacitance measurements (VSCAP) on bne-MeV electron irradiated Si diodes . 
Detailed descriptions on how these parameters were obtained are given in section 
2.5 (page 46) where they are used in the computer calculations. 
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The coupled gold donor-acceptor recombination center was also used in 
TR-1. For high efficiency silicon solar cells and silicon devices with long 
recombination lifetimes, the residue recombination center is likely to be gold 
due to gold contamination during the device processing steps. The two-level 
gold recombination, center is also a good model for other possible residue metallic 
impurities in terrestrial silicon solar cells made from low-cost silicon. 

Experimental evidence from phosphorus diffusion profile studies has 
indicated that the phosphorus-vacancy center is necessary to explain the 

2 

non-complementary error function or non-Gaussian phosphorus diffusion pi^'ofiles. 

A similar defect or defect-impurity cluster with a donor charge state may be 
present in the study of the boron diffusion profile in silicon.^ Although 
these diffusion profile studies and models by no means indicate the presence 
of these two centers, they will be employed to determine if recombination at 
these centers will limit the ultimate efficiency of silicon solar cells when 
all other impurity recombination centers are eliminated. 

In order to simulate the experimental high efficiency silicon_solar cells - 

produced by Sandia Corporation,^ the concentration of the recombination centers 

in the quasi-neutral base region is selected so that the open circuit .midbase 

lifetime is equal to the reported value, 700 ps, in the high-efficiency Sandia 

P+/N/N+ cells. For the gold model, this lifetime corresponds to a gold concen- 

10 3 

tration of 4.465x10 Au/cm , an extremely low value. 

The interband Auger recombination mechanism is also included in the 
computer calculations since it was thought that Auger recombination may be 
important. The Interband Auger recombination rates are from experimental data.-* 

The impact ionization rates are computed from mass action law e^c’^? and 
e^=c^n? where n^(297)=l. 036x10^*^. The interband Auger capture rates and tbeir 
inverse or impact generation rates are given in the following table. 
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Table 2.4 Interband Auger Recombination 
and Impact" Generation Rate of 
Electron-Hole Pairs by Electrons 
and Holes in Silicon at 29 7K. 

Lifetime at 10 carrier /cm 


n 

c 

n 

e 

-31 

= 2.8 X 10 

=3.0 X lO"^^ 

cm^/sec 
1 /sec 

T = 2.6 

P 

ys 

cP 

-32 

= 9.9 X 10 

cm^/sec 

T = 10.1 

n 

US 

eP 

= 1.06 X 10“^^ 

1 /sec 


XI, 

•1 

= 1.036 X 10^° 

1 /cm^ 





Other complicating factors are not included in the computer model, 
following the reasons given in TR-1.^ There are two considerations, one concerns 
the fundamental material physics and the other the availability of reliable data. 

Energy gap narrowing and Fermi-Dirac statistics are both excluded in the 
model for a number of fundamental reasons given below, (i) Energy gap narrowing 
has been attributed to high impurity density which gives rise to random impurity 
potentials. These random potentials perturb the band states into the band gap 
to form the so-called band tail and gap states. However, these band tail and gap 
states are highly localized at the impurity sites. They are recombination or 
trapping states with nearly zero or very low mobilities and hence cannot be 
treated in the same way as the high mobility band states. Thus, they do not 
contribute to energy gap narrowing in the way that is taken into account in the 
literature where an effective intrinsic carrier density, n^, is used which 
depends on an effective energy gap that includes narrowing. Such a treatment 
tacitly assumes that band tail states are high mobility states but they are not. 
(ii) Localized statistics for bound state energy levels, instead of Fermi-Dirac 
statistics for band states, must be used to determine the electron or hole 
occupation at these localized band tail states, (iii) Deionization of the shallow 
energy level impurities (phosphorus, arsenic and boron) must be taken into 
account at high carrier concentrations. (iv) Impurity banding must also be taken 
into account at high impurity concentrations when inter- impurity distance becomes 
comparable to the bound state orbital diameter. 

Although energy gap narrowing, Fermi-Dirac statistics, impurity deioniza- 
tion and impurity band formation are excluded in the computer model, their effects 
on solar cell performance tend to compensate each other, resulting in a much 
smaller net contribution. This compensation effect also makes the experimental 
separation of their influences very difficult. 
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To illustrate the compensation effect, we note that energy gap narrowing 
will increase the minority carrier density and hence increase the recombination 
current. However, Fermi-Dirac statistics will decrease the minority carrier 
density and hence decrease the minority carrier recombination current. 

As another illustration, impurity deionization and banding will decrease 
the majority carrier concentration and hence increase the minority carrier 
concentration, recombination rate and current. However., reduction of the majority 
carrier concentration will reduce the Debye screening of the electron-hole 
Coulombic interaction. Debye screening of the electron-hole attraction is the 
cause of real energy gap narrowing since it involves no bound states and only 
band electrons and holes. Thus, smaller Debye screening due to lower majority 
carrier density from impurity deionization would lessen energy gap narrowing and 
hence decrease the recombination rate and current. 

Aside from the fundamental considerations given above which indicate that 
the net effect of the four factors we neglected (energy gap narrowing, Fermi- 
Dirac statistics, impurity deionization and impurity banding) is small, there 
are also several experimental indications which are not inconsistent with the 
contention that their net effect is small. First, there is no reliable and 
consensus experimental data on energy gap narrowing. Some recently reported 
data were subsequently proved to be inaccurate. Second,, silicon solar cell 
performance data (efficiency, VOC, JSC)' of the best high-efficiency cells are 
already nearly equal to the theoretically computed maximum values without the 
inclusion of energy gap narrowing, Fermi-Dirac statistics, impurity deionization 
and impurity banding. The lack of reliable experimental data suggests that the 
factors are compensating so that it is difficult to isolate them out individually 
from experimental data. The good agreement between theory and experiments on 
the best high efficiency silicon solar cells again suggests that these factors 
are compensating or individually unimportant. 
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2.1.2 Computer Model 


Numerical computations are made using the transmission line equivalent 
circuit model method described in section 2.1 of Technical Report No. 1. The 
circuit model was synthesized from the exact Shockley Equations. Hence, the 
results are exact with only numerical inaccuracy which is much less than 1%. 

The only changes made are in the computer model and the numerical integration 
scheme for the photocurrent generators. 

The computer model and the numerical integration scheme for the photo- 
current employed in this chapter give slightly better accuracy than those 
employed in T.R. No. 1. The computer model is again a four region model 
consisting of the emitter region, the junction space charge .region, the quasi- 
neutral base and the back surface field region. A two hundred (198) section 
model is used in which the section numbers of the four regions are 25, 50, 73 
and 25 in contrast to 20/60/91/25 used in TR-1. The emitter region thickness 

_4 

is taken to equal the emitter junction depth, 0.25x10 cm. The emitter space 

-4 

charge region is chosen to be 0. 75x10 cm which insures that it wil l c over 
the entire transition region of rapidly changing electrostatic potential at 
equilibrixim (zero applied voltage in the dark). The *true* space charge region 
will always be thinner than this second region of the 4-region model. The fourth 
region is made large enough to cover the entire back surface field layer. The 

thickness of the fourth region is set to eight times the emitter junction depth, 

-4 -4 

2.0x10 cm. The physical BSF layer has a thickness of 1.0x10 cm. The quasi- 

-4 

neutral base or the third region has a thickness of 300-0.25-0.75-2.0=297x10 cm. 
The cell thickness is. chosen to be 300x10 ^ cm to make it equal to that of the 
experimental high efficiency cell. The section sizes are constant except in the 
third or quasi-neutral base region where the size increases linearly to a peak 
and then decreases linearly so that the sizes of the boundary sections are 
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identical to those of the adjacent regions, the space charge region and the 
back surface field region. Such a section size variation without disconti- 
nuities across the region boundaries is crucial to obtain accurate numerical 
solutions. 

The numerical integration method for the photocurrent generators in 

each section has been improved over that used in TR-1. In TR-1, the Trapezoid 

rule was used for the double integral of integration over the position in the 

cell and over the wave length of the solar spectra. For the 200 ym thick cell, 

this gave 34.405 mA/cm maximum photocurrent for the 22/60/91/25 computer 

model in TR-1. The exact integration over position (equivalent to ca ! ca j ca ! Qo 

computer model) and Simpson rule integration over solar spectrum used In 

2 

this report gave 33.880 mA/cm . Thus, the TR-1 computer model corresponded 

to a slightly higher solar intensity than that used in the input spectra due 

to the numerical error in the Trapezoid integration rule used in TR-1. We 

2 

shall call 33.880 mA/cm the exact value (i.e., «>/oo/co/oo model) with the 
understanding that it is obtained by a Simpson Rule integration over wavelength 
which makes it still approximate but is dictated by the discrete solar 
spectral irradiance data. 

In the new scheme, the distributed photocurrent generator in each 
section is evaluated by the Simpson rule over the wavelength integration and 
exactly over the position integration. The photocurrent generator located 
at position x- is then given by 

qgo(xi)[x2"X^.I = q| “ G^(0,X)[e~^"'l - (2.1) 

The sum of the distributed photocurrent generators over all 198 sections 

2 

is equal to the exact value of 33.880 mA/cm indicated above. This removes 
an absolute n^lmerical error associated with the use of the Trapezoid rule for 
the dx integration used in TR-1. 
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The test of the two methods of integration on solar cell parameters 

is illustrated in Table 2.5. The cell is a silicon P+/N/N+ BSF cell Td.th 

Auger recombination and spatially constant gold recombination density of 
10 3 

4.465x10 Au/cm to give an open circuit midbase steady-state lifetime of 

Mss^sured phosphorus and boron impurity profiles were 
also used (see next subsection 2.1,3) which give 

^AA^^^ " == 1.25xl0^°[l - (Z/L^)2/3]U(L^_Z) 

- 5.0x10^^ 

- 2.5xlO^®expt-(L-Z)®/L2l Atoms/cm^ (2.2) 
Here, U(X) is the unit step function: 

U(L^— 2)=! for Z<Lj^ and U(L^— Z)=0 for Z>L^. is evaluated to give 
a emitter junction depth of 0.25yia. is evaluated to give a BSF layer 
junction of 1.0 ym from the back surface. 


Table 2.5 Effect of Photogenerator on Silicon P+/N/N+ Cell 
Performance at one AMI Sun 


Computer 
Run No. 

JPHOTO 

mA/cm^ 

EFF 

% 

VOC 

mV 

JSC 

mA/cm^ 

FF 

Integration 

Method 

dX 

Integratio] 

Method 

dx 

RUN571 

35.030 

19.527 

669.40 

33.877 

0,76566 

Trapezoid 

Trapezoid 

RUN575 

34.448 

19.197 

668.68 

33.345 

0.76557 

Simpson 

Exact 


It is evident that the errors are insignificant compared with achievable 
experimental accuracy. However, when comparing various fundamental mechanisms 
using the computer results, the more accurate model given by RUN575 using the 
Simpson integration rule for dX and exact integration rule for dx is preferred, 
CPU time was increased by much less than 1%. 
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2.1.3 Device and Material Parameters 


The following parameters are used so that the results of this study 


on the various recomhination effects can also be compared with published 


experimental data. The experimental parameters of the Sandia P+/N/N+ cells 


are used here. They are listed below. 

Table 2.6 

Item 

1 T = 297 K 

2 n^ = 0.10360 E+11 1/cm^ . 

3 L = 300.00 pm 

4 Xj = 0.25 ym 

5 = 1.0 pm 

6 Pt> “ 10 (n— type FZ) fi-cm 

® 3 

7 Cq = 0.125 E+21 B/cm-^ 

8 = 0.500 E+15 P./cm^ 

^ 3 

9 C= 0.250 E+21 P/cm or 

0.700 E+21 P/cm-^ 

10 N^(x) = Cq[1-(x/L^)^^^]U(L^-x^) B/cm^ 

11 Np(x) = Cj^expt-CL-x)^/!^] P/cm^ 


Reference 


4 

4 

4 

4 

3.4 
4 

2.4 

3,9 

2,8 


The concentration profiles of boron diffused emitter and the phosphorus 

diffused BSF given in items 10 and 11 of this table come from sheet resistivity 

2 3 8 9 

measurements xd-th anodic etching of thin layers. ’ ’ ’ The formulae given 

in items 10 and 11 are simple empirical fits to the experimental data 

2 3 9 

given by the BTL (Bell Telephone Lab) work. ’ The Westinghouse boron profile 
2/3 

also fits (1-Y ) reasonably well. However, the two step phosphorus profile 

g 

reported by Westinghouse (and IBM) and attributed to P^O^N^ nitride layer 

formation is not used in this chapter but is used in chapter 3 where studies 

2/3 

of Ti-doped Westinghouse cells are made. Both the (1-Y ) boron profile 
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g 

and the exp(-Y ) phosphorus profile are considerably flatter than the ideal 

2 

diffusion profiles, exp(-Y ) and erfc(Y) . 

2 

Two values for the phosphorus surface concentration, , are used as 

20 

listed in item 9. The lower value (2.5x10 ) corresponds to the electrically 

20 

active phosphorus or electron concentration while the higher value (7.0x10 ) 

20 

corresponds to the electrically active phosphorus donors, but (4.5x10 ) of 

20 

these may be neutralized by the PV acceptors to give a net of (2.5x10 ) 

3 

electrons/cm if the PV pairs are electrically active. 

The phosphorus surface concentrations given in item 9 and the boron 

2 3 

surface concentration given in item 7, both from BTL work,’ can also be computed 

from the sheet resistance measurements of Sandia^ and Westinghouse^*^ using a 

10 

resistivity versus concentration curve similar but not identical to that of 

11 19 ”3 

Irwing’s. Differences occur for concentrations greater than 10 cm 

where our mobilities drop to constant values of y^=90 and Pp=45 as given by 

Equations (2.60) and (2.61) of TR-1. 

The empirical fits to th e experi men ta l di ffusion pro files of phosp hor us . 

6 2/3 

exp(-Y ), and boron, (1-Y ^) , will be known as the profiles in this 

report . 
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2.1.4 Solar and Optical Parameters 


The AMI solar radiation spectrum used in this report is identical to 

12 

that given on p.38 and 39 of TR-1 which come from a Table given by Thekaekara. 

As in TR-1, the data for a particular ambient condition is used since a 

detailed tabulation of the AMI solar spectrum is given. The ambient corresponds 

2 

to 20mm H^O and 3.4mm 0^ and the total incident power is 88.92 mW/cm . This 

is somewhat smaller than those used by other workers in AMI efficiency calcula- 

2 13 

tions which ranged upwards to 92.46 mW/cm such as the work by Fossura and 

14 

the Westinghouse experimental quartz iodine simulator calibration by a 
NASA Lewis standard. 

Accurate absorption curves in silicon were obtained from original 
1 

publications at 291K. These are fitted and interpolated at closely spaced 
wavelengths. The silicon surface reflection coefficient, R(X) , is calculated 
from 


4n 

R(X) = 

(1+n) ^ 


(2.3) 


Here, the wave length dependence of the index of refraction, n(X), and extinction 

coefficient, k(X) , are read off the graphs on page B-2 and B-3 of the JPL 

15 

Solar Cell Radiation Handbook of Tada and Carter. Reflection is neglected 
in the study of high efficiency cells which have antireflection coating but 
included in Westinghouse bare surface cells. 

These are tabulated in Table 2.7. The four significant figures 
listed do not imply absolute accuracy but serve as a reference in the 
computer model so that the numerical results will retain their relative acduracy-';, , 
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Table 2.7 AMI solar spectral irradi'ance, silicon absorption 

coefficient and silicon surface reflection coefficient 
as a function of photon wave length. 88'. 92 mW/cm^ 

(Part 1, 2900-6000A) 


X(li) 

. 2900D+00 
. 2950D+00 
. 3000D+00 
. 3050D+00 
. 3100D+00 
. 3150D+00 
. 3200D+00 
. 32S0D+00 
. 3300D+00 
. 3350D+00 
. 3400D+00 
. 3450D+00 
. 3500D+00 
. 3550D+00 
. 3600D+00 
. 3650ri+00 
. 3700D+00 
. 3750D+00 
. 3B00D+00 
. 38S0D+00 
. 3900D+00 
. 39500+00 
. -10000+00 
. 40500+00 
.41000+00 
. 41500+00 
. 42000+00 
. 42500+00 
. 43000+00 
. 43500+00 
. 44000+00 
. 44500+00 
. 45000+00 
45S0D+L0.0 
, 46000+00 
. 46500+00 
. 47000+00 
. 47500+00 
. 48000+00 
. 48500+00 
. 49000+00 
. 49500+00 
. 50000+00 
. 50500+00 
. 51000+00 
. 51500+00 
. 52000+00 
. 52500+00 
. 53000+00 
. 53500+00 
. 54000+00 
. 54500+00 
. 55000+00 
. 55500+00 
. 56000+00 
. 56500+00 
. 57000+00 
. 57500+00 
. 58000+00 
. 58500+00 
. 59000+00 
. 59500+00 
. 60000+00 


E('eV) 

. 42750+01 
. 42030+01 
.41330+01 
. 40650+01 
. 40000+01 
.39360+01 
. 38750+01 
.38150+01 
. 37570+01 
.37010+01 
. 36470+01 
. 35940+01 
. 35430+01 
. 34930+01 
.34440+01 
. 33970+01 
. 33510+01 
33070+01 
. 32630+01 
,32210+01 
. 31790+01 
.31390+01 
. 31000+01 
. 30620+01 
.30240+01 
. 29880+01 
. 29520+01 
.29180+01 
. 28840+01 
. 28510+01 
.28180+01 
. 27860+01 
. 27550+01 
,.27250+0.1 
.26960+01 
26670+01 
. 26380+01 
.26100+01 
, 25830+01 
. 25570+01 
,25310+01 
. 25050+01 
. 24800+01 
, 24550+01 
. 24310+01 
, 24080+01 
. 23850+01 
. 23620+01 
.23400+01 
.23180+01 
. 22960+01 
. 22750+01 
.22540+01 
.22340+01 
,22140+01 
,21950+01 
.21750+01 
.21560+01 
.21380+01 
.21200+01 
.21020+01 
. 20840+01 
.20670+01 


a(cm“^) 
.21190+07 
. 18680+07 
. 16820+07 
. 15470+07 
. 14240+07 
. 13440+07 
. 12700+07 
. 12040+07 
. 11510+07 
. 11210+07 
. 11030+07 
. 10890+07 
. 10660+07 
. 10500+07 
. 10030+07 
. 88380+06 
. 61940+06 
.44190+06 
. 29790+06 
. 20440+06 
. 15320+06 
. 12420+06 
. 95960+05 
, 74530+05 
. 63890+05 
. 59030+05 
. 56190+05 
.53410+05 
. 48400+05 
.41110+05 
. 33990+05 
.29010+05 
. 26750+05 
..2.4S4D.+.05. 
.22380+05 
. 20260+05 
18190+05 
. 16160+05 
. 14930+05 
. 14190+05 
. 13420+05 
. 12650+05 
. 12020+05 
.11430+05 
. 10810+05 
. 10200+05 
. 96910+04 
. 92370+04 
.88330+04 
. 84830+04 
. 81250+04 
.77040+04 
. 72900+04 
. 70100+04 
. 67710+04 
. 65250+04 
. 62550+04 
. 59910+04 
. 57790+04 
. 55960+04 
. 54160+04 
. 52080+04 
.49990+04 


. 00000+00 
. 00000+00 
, 41000+01 
, 11400+02 
. 30500+02 
. 79400+02 
. 20260+03 
. 26950+03 
.33160+03 
, 38340+03 
. 43130+03 
. 44920+03 
. 48050+03 
. 49800+03 
. 51370+03 
. 56130+03 
. 60350+03 
. 60940+03 
. 60800+03 
. 60980+03 
. 62390+03 
. 69120+03 
. 84990+03 
. 99280+03 
. 10740+04 
.11050+04 
. 11040+04 
. 10870+04 
. 10680+04 
. 11000+04 
. 12160+04 
. 13100+04 
. 138BD+04 
...1.435D.+04 
. 14520+04 
, 14510+04 
. 14510+04 
. 14700+04 
. 15030+04 
. 14430+04 
. 14350+04 
, 14540+04 
, 14510+04 
. 14400+04 
. 14170+04 
. 13850+04 
. 13900+04 
. 14100+04 
. 14070+04 
. 13940+04 
, 13720+04 
. 13540+04 
. 13370+04 
, 13360+04 
. 13190+04 
. 13300+04 
. 13380+04 
. 13470+04 
. 13470+04 
. 13470+04 
. 13410+04 
. 13290+04 
. 13200+04 


l-R(X) 
33290-00 
.38810-00 
.38810-00 
. 40650-00 
. 41980-00 
. 42850-00 
. 43660-00 
. 44300-00 
. 44650-00 
. 44700-00 
. 44650-00 
. 44120-00 
. 43300-00 
. 42500-00 
. 41800-00 
, 41360-00 
. 41370-00 
, 42200-00 
. 43890-00 
. 45750-00 
. 47420-00 
. 49000-00 
. 50400-00 
.51500-00 
. 52560-00 
. 53550-00 
. 54590-00 
. 55400-00 
.56140-00 
. 56850-00 
. 57500-00 
. 58060-00 
. 58590-00 
..59.120.-00_ 
. 59580-00 
.59970-00 
. 60350-00 
. 60660-00 
. 60970-00 
. 61250-00 
. 61520-00 
. 61760-00 
. 62000-00 
. 62200-00 
. 6240D-00 
. 62620-00 
.62810-00 
. 63000-00 
.63210-00 
, 63380-00 
, 63520-00 
. 63720-00 
, 63880-00 
. 64000-00 
,64110-00 
, 64260-00 
, 64370-00 
. 64480-00 
. 64580-00 
. 64690-00 
. 64770-00 
.64860-00 
. 64970-00 
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- Table 2. 


7. AMI solar spectral irradiance, silicon absorption’ 

coefficient and silicon surface reflection coefficient 
as a function of photon wave length. 88.92 mW/cm^ 

(Part 2, 6000-12640A) 


A(w) 

. 6000D+00 
. 6050D+00 
. AlOOD+00 
. 6200D+00 
. 6300D+00 
. 6400D+00 
. 6500D+00 
. A600D+00 
. 6700D+00 
. 6800D+00 
. 6900D4'00 
. 6980D+00 
. 7000D+00 
. 710011+00 
. 7200D+00 
. 7280D+00 
. 7300D+00 
. 7400D+00 
. 7500D+00 
. 7620D+00 
. 7700D+00 
. 7800U+00 
. 7900D+00 
. BOOOD+00 
. 80<50D+00 
. 8250D+00 
. 8300D+00 
. 8350D+00 
. 8460D+00 
. 8600D+00 
. 8700D+00 
. 8750D+00 
. 8870D+00 
. 9000D+00 
. 9070D+00 
. 9150D+00 
. 9250D+00 
. 93000+00 
. 94000+00 
. 95000+00 
. 95500+00 
. 96500+00 
. 97500+00 
. 98500+00 
. 10180+01 
. 10820+01 
. 10940+01 
. 10980+01 
. 11010+01 
. 11280+01 
. 11310+01 
. 11370+01 
. 11440+01 
. 11470+01 
. 11780+01 
. 11890+01 
. 11930+01 
. 12220+01 
. 12360+01 
. 12640+01 


E(eV) . 

. 20670+01 
. 20500+01 
. 20330+01 
. 20000+01 
. 19680+01 
. 19370+01 
. 19080+01 
. 18790+01 
. 18510+01 
. 18230+01 
. 17970+01 
. 17760+01 
. 17710+01 
.17460+01 
. 17220+01 
. 17030+01 
. 16990+01 
. 16760+01 
. 16530+01 
. 16270+01 
. 16100+01 
. 15900+01 
. 15700+01 
. 15500+01 
. 15380+01 
. 15030+01 
. 14940+01 
. 14850+01 
. 14660+01 
. 14420+01 
. 14250+01 
. 14170+01 
13980+01 
. 13780+01 
. 13670+01 
. 13550+01 
. 13410+01 
. 13330+01 
. 13190+01 
. 13050+01 
. 12980+01 
. 12850+01 
. 12720+01 
. 12590+01 
.12180+01 
. 11460+01 
. 11330+01 
. 11290+01 
. 11260+01 
. 10990+01 
. 10960+01 
. 10900+01 
. 10840+01 
. 10810+01 
. 10530+01 
. 10430+01 
. 10390+01 
. 10150+01 
. 10030+01 
. 98100+01 


a(cm”^) 
499935+04 
. 47950+04 
. 46270+04 
. 43000+04 
39820+04 
. 36750+04 
. 33770+04 
. 31300+04 
.29060+04 
. 27140+04 
. 25290+04 
. 23850+04 
. 23500+04 
.21790+04 
. 20330+04 
. 19200+04 
. 18920+04 
. 17540+04 
. 16200+04 
.14640+04 
. 13620+04 
. 12590+04 
. 11780+04 
. 11000+04 
. 10440+04 
. 87440+03 
. 83580+03 
. 80000+03 
. 72280+03 
. 63060+03 
. 57090+03 
.54160+03 
.47380+03 
.41320+03 
33130+03 
. 34550+03 
. 30530+03 
. 28660+03 
. 24970+03 
.21360+03 
. 19680+03 
. 16990+03 
. 14350+03 
. 11770+03 
. 52000+02 
. 65680+01 
. 47500+01 
.42100+01 
.38500+01 
. 16000+01 
. 14600+01 
. 11500+01 
. 88000+00 
. 76000+00 
. 62000-01 
. 25000-01 
. 19000-01 
. 40000-02 
. 20000-02 
. 00000-00 


P(W/m^U) 

, 13200+04 
. 13110+04 
, 13080+04 
. 12940+04 
. 12810+04 
. 12720+04 
. 12570+04 
. 12440+04 
. 12270+04 
. 12100+04 
. 11960+04 
. 10100+04 
. 11750+04 
. 11570+04 
. 11350+04 
. 10030+04 
. 11180+04 
. 10950+04 
, 10770+04 
. 79400+03 
. 10390+04 
. 10190+04 
. 10000+04 
, 98120+03 
. 87440+03 
.93160+03 
.92180+03 
.91240+03 
. 47620+03 
. 50640+03 
. 45380+03 
. 44920+03 
. 44860+03 
. 44890+03 
. 45520+03 
. 46150+03. 
. 27900+03 
. 22180+03 
. 31340+03- 
. 29650+03 
. 32110+03 
. 34440+03 
. 57690+03 
. 54460+03 
.61750+03 
.51290+03 
. 46410+03 
. 50370+03 
. 50480+03 
. 13510+03 
. 15220+03 
. 14310+03 
. 19120+03 
. 17450+03 
. 39930+03 
. 40220+03 
. 42400+03 
.39180+03 
. 39080+03 
. 32920+03 


,1-R(X) 

. 64970-00 
. 65030-00 
. 65120-00 
. 65290-00 
. 65430-00 
. 65590-00 
. 65710-00 
. 65870-00 
. 66000-00 
. 66070-00 
. 66200-00 
. 66320-00 
. 66380-00 
. 66480-00 
. 66600-00 
. 66700-00 
. 66720-00 
. 66790-00 
. 66900-00 
. 67010-00 
. 67050-00 
, 67100-00 
. 67260-00 
. 67350-00 
. 67380-00 
. 67460-00 
. 67500-00 
. 67S2D-00 
. 67580-00 
. 67660-00 
.67710-00 
. 67740-00 
. 67800-00 
. 67840-00 
67880-00 
. 67910-00 
. 67950-00 
. 67970-00 
. 68000-00 
. 68050-00 
, 68070-00 
.68100-00 
. 68130-00 
. 68130-00 
. 68270-00 
. 68420-00 
. 68450-00 
. 68470-00 
68480-00 
. 68500-00 
. 68520-00 
. 68530-00 
. 68540-00 
. 68550-00 
. 68600-00 
. 68610-00 
. 68620-00 
. 68660-00 
. 68690-00 
68730-00 
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The dependence of the maximum photocurrent or areal photon, density 
absorbed per second on the sample thickness is a, very useful parameter in 
measuring the performance of an experimental solar cell and in projecting 
probable further improvement of a particular technology. For these reasons, 
the results of the Simpson rule integration -of dX and exact integration of 
dx is given in the following table for silicon at one AMI sun. Reflection 
at the silicon surface is neglected. 

Table 2.8 Maximum Photocurrent as a Function of Silicon 
Solar Cell Thickness 


L (pm) 

JPHOTO (mA/cm^) 

0.1 

3.1318 

0.2 

4.7755 

0.5 

8.1773 

1.0 

11.812 

2.0 

16.219 

5.0 

22.296 

10.0 

26.208 

20.0 ■ 

29.088 

50.0 

31.491 

100.0 

32.772 

-_150.0- 

33.438 

200.0 

33.880 

300.0 

34.453 

400.0 

34.826 

500.0 

35.101 

1000. 0 

35.943 

2000.0 

36.839 

5000 

37.702 

10000 

37.971 

CO 

40.853 


2 2 
The difference of 0.005 mA/cm between the L=300ym result (34.453 mA/cm ) and 

2 

that given in Table 2.5 of RUN575 (34.448 mA/cm’ ) comes from neglecting the 
photon flux in the last dX at the absorption threshold. It should be mentioned 
that results of JPHOTO given in the above table represent the maximum short 
circuit current obtainable if the back surface is perfect reflecting. If the 
back surface is perfect transmitting, i.e. one-pass, then JSC (max) < JPHOTO. 
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2.2 


Results on Back Surface Field 


The effect of back surface field (BSF) on the solar cell performance is 
analyzed in this section. Different diffusion profiles would give different 
amount of back surface fields. Its effect on solar cell performance is 
also analyzed. 

The cell structures used here are the P+/N and the P+/N/N+ types 
whose impurity concentration profiles are given by 

(1) P+/N Cell 

••aA - "dD - Cgexpl-CZ/Lph . (2.4) 

(2) P+/N/N+ Cell-2 

Naa - % " - Cg - Cgexp[-(L-Z)^/L2] (2.5) 

(3) P+/N/N+ Cell-6 

Naa - Kgg = C^expE-CZ/L^)^] - - C^exp [- (L-Z) ^/Lg] (2.6) 

o r\ 1 / 0 

Here, C =1.25x10 , C =5.0x10 and C_=2.5xl0 . The gold recombination model 

U B L ’ 

10 3 

is used with a spatially constant gold concentration of N^^=4. 465x10 Au/cm 

to give a midbase open circuit steady-state lifetime of 700 ps in the P+/N/N+ 

cell. This gold density is determined from the steady-state lifetimes versus 

position of a trial computer run in which the gold concentration was assumed 

to be 1.0x10 Au/cm . The midbase steady-state lifetimes at the open circuit 

condition for the- P+/N/N+ Cell-2 are x =718.996 and x =718.987 ys at the gold 

n p 

10 3 

concentration of 4.4654x10 Au/cm . These values are sufficiently close to 
the reported experimental values, 700 ps, so that a second iteration was not 
made. 

A very Important condition used in this calculation is that both the 
front and the back surfaces are short-circuited. The short-circuit boundary 
condition corresponds to infinite surface recombination rate or -velocity. In 
addition, the surfaces are assumed to be neutral, that is, the interface 
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recombination centers at the two surfaces are assumed to be at. the neutral 
charge condition. This neutrality condition corresponds to a band of high 
density of interface states whose neutral Fermi level coincides with the 
silicon Fermi level at the silicon surface. Subsequent calculations, using 
a more concrete interface recombination center model, have confirmed the 
correctness of the present short circuit boundary condition. These calculations 
also showed that interface recombination had very small effect on the solar 
cell performance. This is an important result in view of the considerable 
current research efforts on the improvement of solar cell performance by 
reduction of interface or surface recombination velocity. Our calculations 
show that interface recombination cannot play a significant role because 
of the large retarding electric field from the diffusion profile of the 
diffused emitter and BSF layers. 

2.2.1 Effect of Back Surface Field. 

The substantial Improvement of the open circuit voltage from the 
presence of a back surface field ,in_a„solar -cell— was- first-recognized" and" 
analyzed in detail by Godlewski, Baraona and Brandhorst (GBB).^^ Their 
understanding and analyses were so complete that there is little need to further 
verify the presence of this effect. However, their analytical solutions 
assumed constant impurity concentrations at the two sides of the high-low 
junction barrier. Such an abrupt change of impurity concentration gives the 
upper limit of improvement one can expect from a back surface field layer. In 
practice, the concentration change in the diffused BSF layer is more spread out 
than the abrupt H/L transition assumed by GBB, If the concentration change is 
very gradual, spreading out over many diffusion lengths, then the improvement 
would not be as large. The calculations of this and the next section will 
show that the improvement is indeed very large and there is little difference 

between the two impurity diffusion profiles assumed. 
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The numerical calculations were made using the Trapezoid integration 
rule for the double integral of the distributed photocurrent generator 
before the improvement was made using the Simpson rule and the exact 
integration. They were not repeated using the improvement since the differences 
are small and the results are adequate for the study of the effects from BSF. 

The solar cell characteristics with and without BSF using the 
Gaussian impurity profiles given by Equations (2.4) and (2.5) are compared in 
the following table. 

Table 2.9 Effect of Back Surface Field on Silicon Solar Cell 
Performance at One AMI Sun and 297K. 


Computer 

Cell 

Impurity 

Profiles 

EFF 

VOC 

JSC 

FF 

Run No. 

Type 

Emitter 

BSF 

% 

mV 

mA/cm^ 


RUN555 

P+/N 

exp (-X^) 

none 

16.236 

554.11 

32.456 

.8028 

RUN335 

P+/N/N+ 

exp (-X^) 

exp (-X^) 

19.587 

672.32 

33.954 

.7629 


The F+/N/N+ BSF cell, given by RUN335, shows much higher EFF and VOC but 

only slightly higher JSC than the P+/N cell, given by RUN555. The improvement 

of VOC is large, amounting to 118 mV or almost 4.6(kT/q). This is 21% over the 

2 

P+/N cell which has no BSF. JSC is increased only by about 1.5 mA/cm or about 

4.6% with the BSF. The fill factor, FF, actually decreased from 0.8028 to 

0.7629 which reflects the presence of high injection level condition in the base 

of the BSF cell given by RUN335. The high injection level condition in the base 

also increases the steady-state base lifetime from 562 ps for the P+/N cell to 

719 ys for the P+/N/N+ cell which accounts for the slight increase of JSC just 

noted. The high injection level condition may be further analyzed by noting 

that VOC and FF are directly related if the JSC vs VOC dependence is known. 

For the low injection condition, Jg^=Jj^exp(qVQ^/kT) , while for the high injection 

condition, J„ “J„e3cp(qV-^/2kT) . Thus, for these ideal diode laws, the fill 
oC il UU 
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factor, FF, can be computed for a given These are tabulated below for 

the two cells of Table 2.9 and compared with the exact values given there. 

Table 2.10 Comparison of Exact and Ideal Theoretical Fill Factors. 


Computer 

Cell 

VOC 

FILL 

FACTORS 

Run No. 

Type 

mV 

Exact 

L.L. 

H.L. 

RUN555 

P+/N 

554.11 

.8028 

.818 

.707 

RUN335 

P+/N/N+ 

672.32 

.7629 

CM 

00 

.742 


The results of the P+/N cell given by RUN555 show that the exact FF (0.8028) 

is closer to the ideal low level FF (0.818) than the ideal high level FF (0.707). 

However, a detailed examination of the position variation of the electron and 

hole concentration inside the P+/N cell at the SC, MP (maximxam power) and OC 

conditions indicates that the midbase values of the minority carrier or hole 

13 13 14 

density in the n-base are 1x10 , 2x10 and 2x10 while the equilibrium 

14 -3 

electron density is 5x10 cm . Thus, high level condition begins to set in 
under the open circuit condition in the P+/N cell at one AMI sun but not at 
fhe* short circuit 6r-max'imum power conditions. 

For the P+/N/H+ cell of Run 335, the exact FF (0.7629) is closer to 
the ideal high level FF (0.742) than the ideal low level FF (0.842). The BSF 
high-low junction has nearly completely shielded the infinite surface recombina- 
tion rate at the back surface. This can also be verified by using the 

16 

analytical solutions obtained by Godlewski, Baraona and Brandhorst in the 
limit of infinite back surface recombination velocity and much larger areal 
donor density in the N+ BSF layer than the N- base layer. In this case, the 
diode current -VO It age characteristics can be approximated by 

(2.7) 

10 “3 2 

Using W =300 ym, n.=i0 cm , t =700 us and J =34 mA/cm , then the computed 
SI h oC 
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open circuit voltage Is 

= (2kT/q)log^[Jg^T^/qn^W^] = (2kT/q) log^(5xl0^) = 678 mV. (2.8) 
This compares favorably with the exact numerical value of 672.32 mV. The 
good agreement shows that indeed the low— high (L/H) potential barrier of the 
N/N+ BSF junction is essentially a perfect reflecting barrier which shields the 
back surface recombination loss completely. 
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2.2.2 Effect of Diffusion Profile 


It was just demons tratied that hack surface recombination is almost' 
completely shielded by the presence of the BSE of a Gaussian diffusion ,pro_file. 
Thus, a more abrupt diffusion profile with higher BSE would not be expected to 
further improve the shielding effect of the L/H or N/N+ potential barrier. A 
comparison of the analytical solution of (2.8) with the exact solution of 
Table 2.10 shows that only about 678-672.32=5.7 mV additional improvement can 
be expected from a larger potential barrier or steeper impurity gradient in the 
BSE layer. This anticipation is demonstrated by comparing the results of the 

Gaussian profile in Tables 2.9 and 2.10 with that of the measured phosphorus 

6 2 
profile, exp(-X ), which is steeper than the Gaussian profile, exp(-X ). The 

net impurity concentration for the steeper BSF cell is given by Equation (2.6). 

The computed characteristics of the solar cells with these two BSF profiles are 

compared in the table given below. 


Table 2.11 Effect of Diffusion Profile of the Back Surface Field 
Layer on Silicon Solar Cell Performance at One AMI Sun. 


Computer 

Cell 

Impurity 

Profiles 

EFF 

VOC 

JSC 

FF 

Run No. 

Type 

Emitter 

BSF 

% 

mV 

mA/cm^ 


RUN335 

P+/N/H+ 

exp (-X^) 

exp (-X^) 


672.32 

33.954 

.7629 

RUN568 

P+/N/N+ 

exp (-X^) 

exp (-X^) 

19.592 

673.64 

33.954 

.7616 


The results are as expected. is not changed and is increased 

by 1.32 mV in the more abrupt measured BSF diffusion profile. The small increase 
of V indicates that the ideal Gaussian profile gives about the same solar cell 
performance as the measured phosphorus diffusion profile. The same conclusion 
may also be made for the measured boron diffusion profile which is in between 
the ideal Gaussian profile and the measured phosphorus profile. 
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2.3 Results on Interband Auger Recombination 

Interband Auger recombination of electrons and holes can be an important 
recombination mechanism in regions of high carrier concentration because its 
steady-state recombination rate is proportional to the square of the majority 
carrier concentration and the steady-state lifetime is inversely proportional 
to the square of the majority carrier concentration. They are given by 


R, = c^N^P + c^P^N 

Auger 


(2.9) 

c"n^P 

(n-Sl) 

(2.9n) 

= c^P^N 

(p-Si) 

(2.9p) 

”^Auger " + 2(c^+c^)PN + c^P^] 


(2.10) 

= 1/c^^ 

(n-Si) 

(2.10n) 

= 1/c^P^ 

(p-Si) 

(2.10p) 


Here, c^ is the interband Auger recombination rate of electron-hole pairs in 

which the energy is carried away by an electron and c^ is the rate in which 

the energy is carried away by a hole. The steady-state recombination rate and 

lifetime given above reduce to the simplified forms shown in Equations (2.9n), 

(2.9p), (2.10n) and (2.10p) in extrinsic material when the majority carrier 

concentration is much higher than the minority carrier concentration. For 

example, in p-Si, the electron lifetime [Equation (2.10p)] could be dominated 

by the interband Auger recombination mechanism if the hole concentration, P, 

is very high or the resistivity is very low. Values in Table 2.4 show that 
X8 3 

at 10 carriers/cm , t (n— Si)=2.6 ps and t (p-Si)-lO ys. Thus, interband 

■p n 

Auger recombination in the base can -be the mechanism limiting the ultimate 
performance of silicon solar cell in low resistivity cells if all other 
recombination losses are eliminated. 

It has also been thought that interband Auger recombination is important 
inside the diffused emitter and BSF layers as well as on their surfaces because 
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of the very high majority carrier concentrations in these layers. This is not 

20 

the case unless the majority carrier concentration is greater than about 10 

3 19 

carriers/cra over the entire layer because the diffusion lengths at 10 

3 

carriers/cra (L -2ym and L =6ym) are much larger than the physical thickness 
P n 

of these layers (about 0.25 to 1 ym) . However, interband Auger recombination 
could be important at the surface of the emitter and the BSF layers where the 
carrier concentration is greater than 10 carriers/cm , provided the surface 
layer is sufficiently thick to give high enough Auger recombination current. 
Whether Auger recombination at the surface is important also depends on the 
magnitude of recombination at the interface or surface states. In our computer 
model, we have assumed an infinite interface recombination rate or velocity. 
Thus, the additional Auger recombination at the surface could not affect the 
results . 

Because the interband Auger recombination rate is strongly dependent on 
the majority carrier concentration, we would expect it to be more important 
_in Uie d iffused ^mitt errand. BSE layers - rf -the impurity- diffusion profile's are 
very flat with high concentrations and drop abruptly at the junctions. The 
measured diffusion profiles of phosphorus and boron are indeed more abrupt than 
the ideal or standard Gaussian profile. Thus, it is useful to determine if 
such an abrupt and high concentration profile can significantly increase the 
interband Auger recombination current in the emitter and the BSF layers. For 
this study, we have varied the diffusion profiles to include the more abrupt 
measured diffusion profiles. The results are described below which confirm the 
initial expectation that the interband Auger recombination mechanism in the 
emitter and the BSF layers are not important, even if the interband Auger 
capture rates, given in Table 2.4, are in error and too low by as much as 100. 
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2.3.1 Emitter Layer 


In order to separate the effect of Auger recombination in the emitter 
from that in the BSF layer, we shall first describe the results on cells with 
no BSF layer. P+/N cells are used in this study but the conclusion is equally 
applicable to the 1^+/P cells since the measured phosphorus and boron diffusion 
profiles have similar abruptness and the Auger recombination rates of electrons 
and holes (Table 2.4) are similar in magnitude. The two p-type emitter 
diffusion profiles used are the Gaussian and the measured profiles given by 


^AA- 

”dD “ C„^xp[-(Z/Lph . 

(2.11) 

^AA- 

"dD - - C3 

(2.12) 


where C=l. 25x10^*^, C =5.0x10^^, U(z) is the unit step function and L- is 
U D J* 

determined by the junction depth of 0.25 pm. The gold recombination model is 

10 3 

again used with a spatially constant concentration of N^^=4, 465x10 Au/cm to 

give approximately 700ps base lifetime. The Auger recombination rates used are 
the latest and supposedly most reliable data which are given in Table 2.4. 

Again the Trapezoid integration rule was used in these calculations and 
the results would not be significantly different if the more accurate Simpson 
rule for dA and exact for dx are used in the double integral of the distributed 
pho to generators . Four computer runs are made for the two diffusion profiles 
with and without interband Auger recombination. The results are given below. 


Table 2.12 Effect of Interband Auger Recombination in the Emitter on 
Silicon Solar Cell Performance at One AMI Sun 


Computer 
Run No. 

Cell 

Type 

Emitter 

Profile 

Auger 
Re comb . 

EFF 

% 

VOC 

mV 

JSC 

mA/cm^ 

FF 

RUN555 

P+/N 

exp (-X^) 

NO 

16.236 

554.11 

32.456 

.80275 

RUN556 

P+/N 

exp (-X^) 

YES 

16.226 

554.09 

32.438 

.80275 

RUN566 

P+/N 


NO 

16.238 

554.11 

32.458 

.80280 

RUN567 

P+/N 

j^_jj2/3 

YES 

16.198 

554,03 

32.382 

.80279 
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As_ expected, interband Auger recombination in the diffused emitter has 

a larger effect in the measured emitter profile than the Gaussian profile. 

(i) It reduces the open circuit voltage by 0.02 mV in the Gaussian profile 

(RUN555 and RUN556) and 0.08 mV in the measured boron profile (RUN 566 and 

RUN567) out of 554.1 mV. (ii) It reduces the short circuit current more than 

2 

the open circuit voltage, but still only by a small amount: 0.018 mA/cm in 

2 

the Gaussian profile and 0.076 mA/cm in the measured boron profile out of 
2 

32.4 mA/cm . Thus, interband Auger recombination in the diffused emitter has 
negligible effect. 

2.3.2 Back Surface Field Layer 

Interband Auger recombination in the BSF layer will reduce the open 
circuit voltage. The L/H junction barrier is less effective in shielding out 
Auger recombination, which is closer to the L/H junction position, than the 
interface recombination at the back surface which is further away from the 
junction. However, numerical results tabulated below again shows that the 
reduction of VOC due to interband Auger_ recombination^ in ^the BSF layer .is .ver.y 
small even for the measured phosphorus diffusion which is more abrupt than 
the ideal Gaussian profile. Six computer runs are given below. The results 


Table 2.13 Auger Recombination in the Back Surface Field Layer 
on Silicon Solar Cell Performance at One AMI Sun 


Computer 

Impurity 

Profiles 

Auger 

EFF 

VOC 

JSC 

FF 

Run No. 

Emitter 

BSF 

Recomb 

% 

mV 

mA/ cm^ 


RUN568 

exp (-X^) 

exp(-X^) 

NO 

19.592 

673.64 

33.954 

.76164 

RUN570 


exp(-X^) 

YES 

19.565 

669.72 

33.935 

.76548 

RUN569 



NO 

19,592 

673.64 

33.954 

.76164 

RUN571 

l_X^/3 


YES 

19.527 

669.40 

33.877 

.76566 

RUN335 

exp(-X^)- 

exp(-X^) 

NO 

19.587 


33.954 

.76294 

RUNS 00 



YES 

19.565 

669.62 

33.936 

.76557 


- 30 - 



















































show that interband Auger recombination in the BSF layer decreases the open 
circuit voltage by about 4 mV for the measured phosphorus impurity profile and 
2.7 mV for the Gaussian profile. It also shows that there is very small difference 
in Vq^ between the Gaussian and the measured phosphorus diffusion profiles. 

The AMI efficiency is changed by less than 0.06%. 

2.4 Results on Surface Recombination 

The exact numerical solution just described in section 2.3 showed that 
interband Auger recombination in the emitter and BSF layers has negligble 
effect on the AMI performance of silicon solar cells. Since the results are 
based on exact numerical solution of a model with infinite surface or interface 
recombination velocity at the interface states, the influence of surface 
recombination was not explored. Although exact numerical solutions have been 
obtained which included finite interface recombination rate, the multitude of 
parameters make the exact solutions complicated to analyze. To make the results 
more transparant physically, a simple model will be used for which simple 
analytical solutions can be obtained. This simple model introduces an Important 
physical parameter, the effective recombination velocity of a layer. It is the 
recombination velocity which will produce the same current as that at the 
recombination center in the bulk of the layer. This model will be used to 
assess the Importance of surface recombination at the Interface states on the 
silicon surfaces of a solar cell. 

2.4,1 Emitter Surface Recombination 

The simple model to be used has a constant majority carrier concentration 
in the emitter and hence constant Auger recombination lifetime. Interface 
recombination is characterized by a constant surface recombination velocity, S^. 

The boundary condition for a P+ emitter is then given by 

Jjj(x=0) = qS^[N(0) - Ng] = - qDp[dN(x)/dx]^^Q (2.13) 
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This surface boundary condition at the emitter surface can be used to get 
a solution of the minoirty carrier (electron) diffusion equation in the p-type 
emitter: 


dX,/dx = qD d^N/dx^ = R. = c^P^N = N/t. (2.14) 

N ^ n Auger Auger 

where the low level condition is assumed since the emitter has such a high 
majority carrier or hole concentration, P. The solution is given by the 
ideal Shockley diode equation with a modified saturation current; 


J = jQ[exp(qV/kT) - 1] 

where the saturation current, J^, is given by 


'^0 = 


(2.15) 


(2.16) 


and S_ is the effective recombination velocity defined by 

£j 


(L /D )S + tanh(W„/L ) 

S_ = (DJ1„) ^ ^ 5-^— 

1 + (Lj,/Dj)Sjtanh(¥j/Lg) 


(2.17) 


D is the minority carrier (electron) diffusivity, D , in the emitter. L Is 
£i n E 

the minority carrier (electron) diffusion length, /d t , where t is the Auger 

n n n 

lifetime given by Equation (2.10^) ._W^_i^t^ Uiickness of t he emitt.er_J.ayer... ~ — 
The saturation current, expressed in terms of the effective recombination 
velocity given above, is particularly useful to determine the effect of interface 
recombination. A sketch of the solutions of this effective recombination 
velocity is given below where S„ is normalized to the diffusion velocity, 

^DF^^E^'^E’ ^E ^E^^DF ~ ^E^ ^^E^^E^ ' (2.18) 

There are four characteristic velocities in this general problem; 


'^Rco 

(ii) v^Q = De^E^^E = ^E^'^E 


(iii) Vj^p ^E^'^E 

(iv) St • 


Recombination Velocity, thick emitter 
Recombination Velocity, thin emitter 
Diffusion Velocity, thin emitter 
Interface Recombination Velocity 


(Wj,«Lj.) C2.19) 
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0 

log I = log(Lg/Wg) 


Figure 2 . 1 The effective recombination velocity of emitter 

recombination current as a function of minority carrier 
diffusion length with surface recombination velocity 
as the constant parameter. 

The three asymptotic lines in the above figure, labled by (1), (ii) and (iii) , 
correspond to the three characteristic velocities. When the emitter is thin 
compared with the minority carrier diffusion length, the effective recombination 
velocity will depend on the surface recombination velocity at the interface. 

This is in the shaded region in the above figure and it is given by 

Sg = (Si + L'h/Q + Sj) 

Sj, - [Sj + (L^/DjWjj)]/[l + <Sj.Wj,/Ilg)]. (2.20) 

Here, Sx-St/v_ =8^./ (D„/W_) . When the emitter becomes very thin, » 

J_ J. Di! X if* i!i th su 

the effective recombination velocity approaches 

= Sj/Cl + Sj) 

Sg = S^/[l + (SjW^/Dg)]. (2.21) 
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Numerical values are computed for the three characteristic velocities 
as well as the corresponding saturation currents for both n-type and p-type 
emitter of a thickness of W =0.25ym. Three constant majority carrier 
concentrations in the emitter are assumed. The numerical results are given 
in Table 2.14. 

Consider the results for thin emitters which are items labeled (ii) 
and (iii) in the table, corresponding to the two asymptotic lines with the 
same labels in the figure on the preceding page. It is evident that the 
effective recombination velocity and the saturation current for the case of 
zero interface recombination, 5^=0, labeled (ii) are very small, giving 
unrealistically high open circuit voltages. 

As the surface recombination velocity at the interface increases, the 
saturation current also Increases and the open circuit voltage decreases. For 
the limiting case of infinite interface recombination velocity, the open 

circuit voltages are given by the line labeled (iii) in Table 2.14. The value 

18 3 19 3 

for the p-type emitter (585 mV at 10 holes/cm and 666 mV at 10 holes/cm ) 

can be compared with the experimental value of the high efficiency Sandia 

P+/N/N+ cell^. Since the actual cell does not have a constant hole concentration 

as we have assumed in the simple model, a comparison is made based on the total 

13 2 

hole concentration in the emitter, Q which is about 20x10 q/cm given in the 

ir"r 

18 3 19 ' 3 

table. This is midway between our two cases of 10 holes/cm and 10 holes/cm . 

A t theoretical V^^ can be computed for the Sandia cell by using an average hole 

18 

concentration, computed from the experimental Qp_j_, which is P^^^=Qp_^/qWg=8xl0 

3 2 

holes/cm . The electron diffusivity at this acceptor concentration is 3.5cm /s 

which gives a diffusion velocity of v =D /W =140Kcm/s. The saturation current 

Dr Tl E 

-13 2 

is 2.8x10 A/ cm and the computed open circuit voltage is then V_ =(kT/q)x 

—13 

log (0.03/ 2. 8x10 )=656mV. Experimental value is 620mV and the difference is 

attributed to recombination in the base layer and in the BSF layer since 

-34- 



Table 2.14 Effects of Surface Recombination on the Interband 
Auger Recombination Current in Thin Emitter 


3 

Pp or (Carriers/ cm ) 

lO^S 

19 

10 

O 

CM 

O 

H 

Wp (}jm) 

0,25 

0.25 

0.25 

Qp (10^\/cm^) 

2,5 

25 

250 

2 

D (cm /sec) 

n 

7.1 

3.0 

2.2 

T (sec) 

n 

10 

100ns 

Ins 

L (pm) 

n 

90 

5.8 

0.5 

(i) (cm/sec) 
(ii) VpQ (cm/sec) 
(iii) v^p (cm/ sec) 

800 

2.5 

284K 

5K 

250 

120K 

44K 

25K 

88K 

(i) (10"'^‘*A/cm^) (Thick) 

(ii) JpQ (Sp=0) (10~^Vcm^) (Thin) 

(iii) Jpp (Sp=<») (10 ^‘^A/cm^) (Thin) 

1.3 

0.0035 

450 

0.8 

0.035 

19 

0.7 

0.35 

1.4 

(i) (mV) (Thick) 

(ii) Vq^ (Sp-0) (mV) (Thin) 

(iii) V_„ (S =■») (mV) (Thin) 

UL» i. 

736 

889 

585 

748 

829 

666 

752 

770 

734 

Qp^(Exp) (lO^^q/cm^) 


20 


V^^(Exp) (mV) (P+/N/N+) 


620 


2 

D (cm /sec) 

P 

2.4 

1.4 

■■■ 

T (sec) 

P. 

2.6ps 

26ns 

0.26ns 

bp (pm) 

25 

1.9 

0.17 

(i) v_ (cm/sec) 

K.00 

(ii) Vpg (cm/sec) 

(iii) Vpp (cm/sec) 

960 

9.6 

96K 

7.4k 

960 

56K 

65K 
■ 96K 
44K 

(i) (10~-'-'*A/cm^) (Thick) 

(ii) (S =0) (10“l'+A/cm2) (Thin) 

(iii) Jjjp (Sp=») (10~l'»A/cm ) (Thin) 

1.5 

0.015 

150 


1.0 

1.5 

0.7 

(i) Vq^ (mV) (Thick) 

(ii) Vq^ (Sp=0) (mV) (Thin) 

(iii) Vq^, (Sp=») (mV) (Thin) 

732 

851 

613 

738 

792 

686 

742 

732 

752 




















































maximum emitter recombination has been assumed in this sample calculation. 

Another useful result from Table 2.14 is that it provides an estimate 

of the range of interface recombination velocity within which significant 

19 

effect can be expected. For example, consider the 10 p-type emitter whose 

data are given in the upper half of Table 2.14. The diffusion velocity, 

line (iii) , is 120,000 cm/sec. Thus, the interface recombination velocity, 

Sj, must be significantly smaller than this value to reduce the recombination 

current. On the other hand, there is no use to reduce the interface recombination 

velcoity to a value below 250 cm/sec, which is the effective recombination 

velocity due to interband Auger recombination in the emitter, 

The simple model allows us to make a simple and rapid estimate of the 

importance of interband Auger recombination in the emitter and at what 

interface recombination velocity does interface recombination become important. 

The parameters that determine the relative importance are the minority carrier 

diffusion velocity, v^^^, = 3nd the effective interband Auger 

2 

recombination velocity in the entire emit.ter., v,,„. = _D„W„/-L -= -W=/t: . -This- • - ' 

effective recombination velocity in the entire emitter not only applies to 

the Auger recombination mechanisms but also the Shockley-Read-Hall thermal 

recombination mec hanism as well .as other localized recombination mechanisms 

as long as a constant recombination lifetime can be defined. As another 

illustration, if the emitter recombination lifetime from all mechanisms is 

1 ns (similar to the third column of Table 2.14), then the effective emitter 

—4 —9 

recombination velocity is v^Q=Wg/T^=0. 25x10 /lO =25K cm/sec. Thus, it would 
be futile to decrease the surface recombination velocity at the interface to 
less than 25K cm/sec. 

These numerical results also show that the conclusion will not be altered 
unless the interband Auger capture rates are in error and too small by 100. 
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2.4.2 Back Surface Field Recombination 

The relative importance of recombination in the bulk of the back surface 

field layer and at the interface states on the back surface cap be estimated 

using the simple model of constant majority carrier concentration employed 

in the emitter in the preceding section. The analytical solutions of Godlewski' 

can be used for the interband Auger recombination mechanism in the BSF layer 

since it has a constant recombination lifetime when we assume -a constant 

majority carrier concentration in this layer. 

The conclusion obtained from the analytical solution is very similar to 

that just obtained for the emitter region:- for a total electron concentration 

15 2 

in the back surface field layer, less than about 10 q/cm and an 

19 

average electron or majority carrier concentration of less than about 10 
3 

carriers/cm , Auger recombination cannot be very important. The relative 
importance of Auger or other recombination mechanisms in the BSF layer to the 
surface recombination at the interface states on the back surface is 
determined in the same way as that in the emitter which was just described in 
section 2.4.3. Thus, the interface recombination velocity on the back surface 
must fall below the diffusion velocity but above the effective recombination 
velocity of the entire BSF layer in order to have any influence on the open 
circuit voltage. 
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2.5 


Results on Enhanced Impurity Solubility 


The, diffused emitter and hack surface field layers are frequently used 
to getter -the recombination impurity centers in the bulk of silicon device? 
because of the enhanced solid solubility in these high concentration layers 
from impurity- shallow donor or impurity-shallow acceptor interactions at the 
diffusion temperatures. In solar cell and other shallow emitter structures, 
the back surface field layer is a particularly effective getter since it is 
formed by high temperature diffusion for a sufficiently long time to allow 
the recombination Impurities in the bulk to outdiffuse-into this layer. 

One common practice is to form the N+ BSF layer at a high temperature to getter 
out all the metallic recombination impurities in the bulk of silicon. 

Gettering by the BSF or diffused emitter layer will improve. the base 
lifetime but at the same time will reduce the emitter and the BSF lifetimes. 

If the emitter or the BSF layer has a high concentration of the recombination 
impurities from gettering action, the short circuit current and the open 
circuit voltage will suffer. The_magnitude_ of this degradation .can he 
predicted by the computer model and is the subject of this section. 

In the following two subsections, we shall describe the results of 
enhanced solubility of the recombination impurity in the emitter of a P+/N 
cell and in the back surface field layer of a P+/N/N+ cell. The two-level 
gold recombination model will be used. The solubility of gold is assumed to 

be proportional to the concentration of the boron or phosphorus impurity in 

18 XO 

the diffused layer 7 A constant background gold concentration of 4.465x10 

3 

Au/cm is also assumed to be present in the entire sample in order to 
simulate the experimental OC base lifetime of 700 ps at one AMI sun reported 

4 

for the Sandia high-efficiency cells. 
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The difference between this study arid that reported in TR-1 is as 
follows. In TR-1, the gold concentration in the base was raised in the same 
proportion as that in the diffused emitter or BSF layer. In this report, 
the gold concentration in the base is kept constant while it is increased 
in the emitter and the BSF layers. The constant gold concentration in the 
base is a better model for the effects of gettering and enhanced solubility 
of recombination impurity centers in the emitter and the BSF layers. The 
gold concentration in both this and previous reports is assumed to be 

proportional to the diffusion impurity concentration. Other possible 

2 3 2 3’ 

dependences are or etc. and or etc. depending 

on the detailed impurity pairing and complexlng mechanisms involved. These 

stronger dependences will decrease the lifetime in the emitter or BSF layer 

near the junction further and hence increase the recombination current even 

more than the linear dependences, and ^AA’ computer 

results of this section show already very substantial reduction of solar cell 

performance from the linear dependences. 


- 39 - 



2.5.1 Emitter Layer 


The principal effect to be expected from an enhanced solubility of 
gold in the diffused emitter is a reduction of the short circuit current 
because of the high absorption of the photons in the emitter layer. Reduction 
of the open circuit voltage should be smaller than that of the short circuit 
current since VOC is mainly determined by base recombination and BSE in this 
high resistivity cell. The first series of study was made to determine the 
effect of emitter impurity diffusion profile since a steeper or more abrupt 
profile would give a higher gold solubility and hence shorter lifetime, closer 
to the emitter junction. This would increase the emitter recombination 
current. The second series of study was made to determine the amount of 
performance degradation as a function of increasing gold solubility in the 
emitter. These are summarized in the following tables. 

The effect of the emitter impurity profile is studied using the 

2 2/3 

standard Gaussian profile, exp(-X ), and the measured boron profile, 1-X ' . 

The cell characteristics are tabulated below. 

Table 2.15 Effect of Emitter Diffusion Profile on Enhanced 
Solubility of Recombination Center and Silicon 
P+/N Solar Cell Performance at One AMI Sun 


Computer 
Run No. 

Emitter 

Profile 

%Au In 
Emitter 

EFF 

% 

VOC 

mV 

JSC 

mA/cm^ 

FF 

RUN556 

exp (-X^) 

0 

16.226 

554.09 

32.438 

.80275 

RUN567 

ragH 

0 

16.198 

554.03 

32.382 

.80279 

RUN564 

exp (-X^) 

1 

15.788 

552.75 

31.643 

.80268 

RUN563 


, 1 

14.951 

551.34 

30.053 

.80237 

RUN565 

exp(-X^) 

10 

15.027 

548.17 

30.362 

.80282 

RUNS 60 . 

1 - 

10 

14.176 

548.50 

28.658 

.80194 
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From the results of the table, we note two important trends. First, 

the enhanced solubility of gold in the emitter can significantly reduce the 

short circuit current but has little effect on the open circuit voltage. 

Second, the diffusion profile also has a significant reduction effect on the 

short circuit current and little effect on the open circuit voltage. The 

2/3 

more abrupt drop-off of the measured boron profile, 1 - X , gives a higher 

gold concentration and lower electron lifetime near the emitter junction than 

2 

the Gaussian profile, exp(-X ). For example, the short circuit current 

2 

decreases from 32.438 mA/cra (RUN556) for the Gaussian emittep impurity 

2 

profile without solubility enhancement to 28.658 mA/cm (RUN560) for the 

measured boron profile with a gold solubility equal to 10% of the boron concen- 

2 

tration. This is a drop of 3.780 mA/cm or nearly 12% in JSC due to enhanced 
solubility of recombination center in the emitter. The drop of JSC accounts 
for the decrease of efficiency from 16.226% to 14.176% since the open circuit 
voltage and the fill- factor both decrease only slightly: AVOC=-5.59mV=-1.0% 
and AFF=-0.00081=-0.10%. 

The results of these examples show that a 2% decrease of AMI efficiency 
may occur, from 16 to 14%, if the recombination center solubility in the 
emitter is enhanced to a value of 10% of the emitter Impurity concentration. 
Such a high value is not unreasonable since impurity interaction can make the 
recombination center equal to dopant concentration under favorable conditions. 
Although the sample calculation is made for a boron diffused emitter of a P+/N 
silicon solar cell, the same kind of performance degradation would be expected 
in a N+/P cell due to a similar behavior of enhanced solubility of the 
recombination centers in the N diffused emitter. 

A more detailed calculation has been made on the performance 
degradation as a function of the gold solubility in the diffused emitter 
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in a P+/N silicon solar cell using the measured boron diffusion profile, 1-X^^^. 

10 3 

A residue gold concentration of 4.564x10 Au/cm is again assumed. ■ Auger ‘ 
recombination is included but its effect is negligible. The results are 
tabulated below. 

Table 2.16 Effect of Enhanced Solubility of Recombination Center 
in the Diffused Emitter on the Performance of Silicon 
P+/N Solar Cells at One AMI Sun 


Computer 

% Au In 

EFF 

VOC 

JSC 

FF 

Run No. 

Emitter 

% 

mV 

mA/cm^ 


RUNS 6 7 

0 

16.198 

554,03 

32.382 

.80279 

RUN559 

.00893 

16.169 

553.98 

32.328 

.80278 

RUNS 6 3 

1.0 

14.951 

551.34 

30.053 

.80237 

RUNS 6 2 

5.0 

14.347 

549.47 

28.946 

.80209 

RUNS 60 

10.0 

14.176 

548.50 

28.658 

.80194 

RUN561 

20.0 

14.030 

547.08 

28.440 

.80181 


An additional result over the previous table is illustrated in this 
„table. It., shows, that, .perjfomance degradation, .tends. ,to slow .down as.„the gold 
solubility increases above the 5 to 10% range. This may be used to simulate 
the operating life of a cell as affected by the diffusion of metallic 
impurities into the emitter from the contact electrodes. The above table 
suggests that there could be a significant initial rapid drop of the JSC due to 
diffusion of the contact metal into the diffused emitter but subsequent changes 
of JSC would be both small and slow as suggested by the saturation trend in 
Table 2.16 wheri gold concentration in the emitter is increased. 

In practice, it has been generally observed that a K+ layer with a high 
concentration of phosphorus is an excellent getter for metallic recombination 
impurities at high temperatures while a similarly strong gettering action has not 
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been observed unambiguously in high concentration boron layers. Thus, the 
results just obtained for the enhanced solubility of gold in the emitter is 
probably more applicable to phosphorus diffused emitter of a N+/P or a 
N+/P/P+ silicon solar cell. < To avoid performance degradation due to 
enhanced solubility of metallic recombination impurities in the diffused 
emitter, one would prefer an arsenic emitter which is probably not as good 
a getter as phosphorus for a N+/P cell or a boron emitter which is likely to 
be a poor getter compared with phosphorus for a P+/N cell. 

The above quantitative result of enhanced solubility of metallic 
impurity in the emitter with enhanced emitter recombination is consistent 
with some experimental observations on the Westinghouse silicon cells and 
on the high efficiency Sandia silicon cells. In the Westinghouse experiments,, 
it has been observed that the Ti-doped P+/N cells have higher efficiency 
than the Ti-doped N+/P cells. Similarly, higher efficiency seems to have 
been obtained in the P+/N/N+ cells than in the N+/P/P+ cells in the Sandia 
high efficiency cell experiments. Base recombination may also be different 
in these cases which may account for the difference between the N+/P and 
the P+/N cells but the results of the emitter recombination enhancement due 
to enhanced solubility in the N**" layer are not contradictory to these 
experimental results. In the Westinghouse cells doped with the highest Ti 
concentration, there is an indication that enhancement of Ti solubility in 
the phosphorus emitter may account for a large part of the efficiency and 
short circuit current degradation. This is described in sections 3,2 and 3.3. 
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2.5.2 Back Surface Field Layer 


Enhanced solubility of recombination impurity in the high - concentration 
diffused back surface field layer would reduce the open circuit voltage because 
of the increased recombination rate and current in the back surface field 
layer. It would have little effect on the short circuit current which is 
mainly affected by recombination near the front surface of the cell. The 
results of computer calculation using the Simpson dl and exact dx integration 
rules for the distributed photocurrent -generators are tabulated in the 
following table. Again, a constant background gold concentration of 
4.465x10 Au/cm is assumed in the entire cell. The gold concentration 
is enhanced in the BSF layer and the enhancement is assumed to be a constant 
percentage of the phosphorus impurity concentration. Auger recombination is 
not included in this computation but its effect was demonstrated to be 
negligibly small in section 2.3. The -measured boron and phosphorus diffusion 
profiles are used in the calculation. .The '-net ionized impurity concentration 
using Table 2.5 is given by 

Nj^(Z) - Njj^(Z) = 1.25xl0^°[l - (Z/Lj^)^''^]U(L^-Z) 

14 

- 5.0x10 

- 2.5xl0^^exp[-(L-Z)^/-L2] lons/cm^ (2.22) 

The results of seven computation runs for increasing gold solubility in the 

BSF layer are tabulated below. They show a significant reduction of the open 

circuit voltage when the gold concentration in the diffused BSF layer is only 

1% of the phosphorus concentration. Tn fact, the values obtained in RTJN576 

4 

are very close to the values of the experimental high efficiency cells. 
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Table 2.17 Effect of Enhanced Solubility of Recombination Center 
in the Back Surface Field Layer on the Silicon F+/N/N+ 
Solar Cell ^performance at One AMI Sun 


Computer 
Run No. 

% Au in the 
BSF Layer , 

EFF 

% 

VOC 

mV 

JSC 

mA/cm^ 

FF 

RUN577 

0.0 

19.261 

672.86 

33.421 

.76162 

RUNS 7 6 

1.0 

18.784 

627.82 

33.373 

.79718 

RUNS 78 

5.0 

18.125 

602.17 

33.242 

.80511 

RUNS 7 9 

10.0 

17.730 

591.01 

33.114 

.80557 

RUN580 

20.0 

17.279 

580.16 

32.913 

.80462 

RUNS 81 

50.0 

16.620. 


32.502 

.80296 

RUN582 

P+/N/«> 

15.991 

553.68 

31.989 

'.80282 


The enhanced solubility of recombination impurities in the BSF layer 
can also be used to model the degradation of operating life due to diffusion 
of metallic impurities from the back contact metal into the BSF layer. From 
the results of RUN577,_which has no gold in the BSF layer, and RUN576, which 
has 1% gold in the BSF layer, it is evident that a large degradation of the 
open circuit voltage can be expected in a short time during the initial 
operation of the cell. Subsequent reduction is considerably smaller for the 
same amount of increase of gold concentration. In the limit when the BSF 
layer is completely saturated with a very high concentration of recombination 
impurities, the performance of the P+/N/N+ cell is degraded to that of a P+/N 

cell (RUN582 on the last line of Table 2.17) if the recombination impurities 
do not penetrate significantly into the base region of the cell. For thick cells, 
such as the 300 pm assumed, deep penetration into the quasineutral base layer is 
probably rather slow compared with the penetration through the 1 ym thick BSF layer. 
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2.6 


Defect-Impurity Recombination Centers. -^ • , 

If all of the metallic Impurity recombination centers can be eliminated, 
then the ultimate efficiency of a silicon solar cell would be determined by 
the remaining deep level recombination cenfers formed by physical lattice 
defects such as silicon vacancies and interstitials; defect complexes such as 
the divacancies and higher vacancy -clusters ; defect- impurity complexes such as 
the oxygen-vacancy pair or A-center, the phosphorus-vacancy pair or the E-center, 
and a boron-vacancy complex yet to be positively identified; shallow-level doping 
or junction forming impurities, such as a -phosphorus or boron pair or clusters; 
and diffusion- induced strain and dislocations. These imperfections can be generated 
in the diffused emitter and BSE layers by the high concentrations of diffusion 
impurities (boron and phosphorus), which form the emitter and BSE layers. 

The computer studies made in this section are aimed to determine the 
silicon solar cell performance limitations imposed by recombination at these 
residual imperfection centers in the diffused emitter and BSE layers. These 
studies are undertaken using two recombination models to simulate recombination 
in the emitter and BSF layers separately. The recombination parameters of these 
recombination models were tabulated in Table 2.2 and 2.3, The basis for selecting 
these parameters are described in the following subsections, 

(A) Emitter Recombination Model 

Eor the diffused emitter layer, the divacancy or boron-vacancy center is 

used as the recombination model. A single level donor-like recombination center 

is assumed whose concentration is given by 

N^,j,(Z) = (FRAC)xl. 25x10^° [1 - (Z/L^)^^^]U(Li-Z) 

+ Centers/cm^ (2.23) 

The factor, FRAG, is the fraction of the boron in the diffused emitter which 

has formed the recombination complex. is the residual recombination density 

1 1j5 
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11 3 

in the entire cell. It covers the range of 0 to 10 Centers/cm and is assumed 

to he spatially constant. The background recombination density can be converted 

into a steady-state base lifetime which can be used as a universal or 

recombination-species-independent parameter under some restricted conditions.^ 

The boron-defect complex has neither been positively identified nor 

characterized. However, transient capacitance measurements on one MeV 

irradiated N+/P silicon diodes have indicated a strong donor recombination 

19 

level at E^+354 mV in the lower half of silicon energy gap. This level was 

r 

also shown to be the main contributor to the leakage current under reverse bias . 
Thus, it will be used as the recombination model in the boron diffused emitter. 
The recombination parameters of this donor level are listed in Table 2.2 and it 
will be referred to as the. boron vacancy complex, although it may involve 

an unknoxTO impurity. 

The recombination coefficients listed in Table 2.2 are obtained as 

follows. first computed from the experimental transient capacitance 

data and e^^ is first computed from the experimental leakage current data. 

Experimental accuracy was two to three significant figures , More significant 

figures are listed to insure that mass action law holds numerically. 

c^Q are estimated from the pre-exponential factor of the thermal emission 

rates, e^^ and e assuming an Arrhenius relationship of the form 
2 

AT exp(— AE/kT). The emission rates are then recomputed so that the mass 

t t tl t 2 

action law at equilibrium .(dark and zero applied voltage), 

is satisfied. The large ratio, ^^-l^^pO ~ consistent with the donor- 

like charge state of this center. Although values, of these parameters given 
in the above table may not be accurate, the solar cell results obtained here 
can also be stated in terms of the electron and hole steady-state lifetimes 


in the base and emitter instead of the concentration of the donor recombination 
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center. Thus, the results of this study on the solar cell performance will not 
be significantly changed when accurate recombination rates become available since 
it is necessary only to add a scale factor to the recombination level density. 

(B) BSP Recombination Model 

The existence of the phosphorus-vacancy or E-center in the diffused 
BSF layer is probably more certain than that of the boron- vacancy complex in 
the diffused emitter. The PV pair has been positively identified by many 
experiments on phosphorus-doped silicon. To study the effect of the PV or 
E center on the silicon solar cell performance, the following position 
variation of the concentration of the PV recombination center is assumed. 


+ (FRAC)x2.5xlO^*^exp[-(L-Z)^/L2l Centers/cm^ (2.24) 

A constant background, denoted by N „ is assumed to simulate the base lifetime 

i. iJ5 

observed experimentally. The factor, FRAC, is the fraction of the phosphorus 
atoms in the diffused BSF layer which have formed the PV pair. Detailed 
-theoretical- analysis- of—experimental "phorsphortTs liif fusion pf o'f lies “have 'sugge¥ted 


a concentration of PV pair with a value of FRAC=0.6 for certain diffusion 
2 

conditions. In this study, FRAC will be varied from 0 to 0.8 to simulate 

a wider range of possible conditions. 

From transient capacitance data, the following recombination parameters 

19 

were obtained for the PV center and are used in this study. The thermal 


emission rate of trapped electrons, e ^i , comes from the experimental data. 
The thermal capture rate, c^q> is computed from the pre-exponential factor 
of the thermal emission rate data using AT exp(-AE/kT). c^^^ is estimated by 

assimiing a hole capture radius of 25A at the negatively charged phosphorus- 

t t 2 “8 2 7 •**6 3 

vacancy pair, giving Cp^=ffpj^0p=Trr^0p=ir(25xlO ) xlO -2x10 cm'^/sec. Such 

a large capture radius is acceptable in view of the large physical size of 
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a pair (it occupies two lattice sites) and the long range Coulomb force 


associated with the negatively charged acceptor state of the PV center. The 


thermal emission rate of trapped hole, e is then computed from the mass 


t 

-pO’ 


action law, e^„ = n^c^„c^, /e^, . The parameters are tabulated in Table 2.3. 
. pO 1 nO pi nl 

Kote that the ratio of the capture rates for the PV center is about 


Cpj^/c^Q = 7.3 compared with the corresponding ratio, c^ l^*^pO ~ 64.3, for the 
boron-vacancy-lmpurity complex listed in Table 2.2. This difference gives 


another variation for checking the sensitivity of the results on the 
recombination model parameters used in the calculation. 
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2.6.1 Emitter Layer 

The effect of recombination at the boron-vacancy center where x may 

be zero) in the emitter layer is more complicated than that at the phosphorus- 
vacancy center in the back surface field layer. This complication arises 
because emitter recombination can affect both J and V „ while recombination 
in the BSF layer reduces mainly and has little effect on 

The reduction of with increasing emitter recombination is large if 
emitter recombination dominates over base recombination. The amount of 
reduction of J with increasing emitter recombination, however, is not 
affected by base recombination, but increasing base recombination will further 
reduce 

A family of curves is computed to show the effect of emitter recombination 

on the efficiency at one AMI sun in silicon P+/N/N+ cells. These curves are 

given in Figure 2.2 where the efficiency is plotted as a function of the 

boron-vacancy complex density, represented by the factor FRAG in Equation (2.23), 

in the_l^ron diffused, .emit ter-.— -FRAG--- T.0‘ 'would"give aT recombination center 

20 3 

density of . 25x10 centers/cm at the front surface. On the x-axis, 

values of two other parameters are also given. These are the open circuit 
minority carrier (hole) lifetime and diffusion length. Since they both vary 
with position, their values at 0.2]jm from the front surface or 0.05ym from 
the P+/N junction are given. 

The density of the base recombination center, simulated by the 


ao 


center, is the constant parameter in these curves, with values of N^^(base)=10 
(curve 10) to 10 (curve 13) centers/cm , Reference curves for negligible base 
recombination are also given which are labeled 1 and lA, Interband Auger recom- 
bination is not included except in curve lA since its effect is small as indicated 
by curve lA. compared with curve 1. The steady-state minority carrier base 
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Figure 2.2 AMI efficiency as a function of boron-vacancy defect 

recombination center concentration and emitter diffusion 
length and lifetime in P+/N/N+ silicon solar cells with 
the base recombination center density as the constant 
parameter. 
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lifetimes are shown next to each curve, ranging from 5.5ys to SSOys. There is 

little variation of the base lifetime at different emitter recombination 

densities along each curve. This is due to the small hole capture rate which 

controls the recombination rate. The much larger electron capture rate 

arises from the Coulomb force of the positively charged donor recombination 

center which gives a ratio of c ,/c „ = 64.3 from Table 2.2. 

n-1 pU 

The effect of emitter recombination on the efficiency can be examined 
by considering the EFF vs FRAC curve of low base recombination, such as curve 
11 which has a base lifetime of 550ps. It is evident that a substantial 
decrease of the efficiency occurs only when the emitter lifetime of electrons 
at the edge of the junction space charge layer becomes substantially smaller 
then 1 ns. A 4% decrease, from 19% to 15%, would occur if the emitter 
recombination center density is 10% (FRAC=0.1) of the boron concentration in the 
emitter. 

For higher base recombination, such as curves 12 and 13 with 55ps and 
5.5ps base lifetimes, the effect of emitter recombination is reduced by about 2 
or a efficiency drop of only about 2%. When base recombination is important, 
increasing emitter recombination will decrease mainly the short circuit current 
which results in efficiency degradation. This is illustrated in the following 
table for the 5.5ys (curve 13) and 550ys (curve 11) base lifetimes. 
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Table 2,18 Effect of the Donor Recombination Center in the Diffused Emitter 
on Silicon P+/N/N+ Solar Cell Performance at One AMI Sun, 


logNTTB= 

11 

13 

11 

13 

11 

13 

11 

13 

FRAG 

EFF(%) 

VOC (mV) 

JSC (mA/ cm^ ) 

FI 

7 

10~^ 

19 -.086 

13.807 

659.58 

516.85 

33.420 

30.042 

.76993 

. 79066 

lO""^ 

18.543 

13.415 

651.97 

516.05 

32.617 

29.239 

,77536 ' 

.79062 

io“^ 

17,184 

12.456 

633.07 

513.94 

30.646 

27.264 

.78759 

.79047 

lO"^ 

16.174 

11.927 

607.30 

512.38 

29.587 

26.203 

. 80043 

.78993 

10“^' 

15; 232 

11:668 

577.17 

510.33 

29.153 

25.767 

.80497 

.78898 

0.8 

13.633 

11.271 

526.30 

498.76 

28.933 

25.547 

.80287 

.78661 


Consider first the short circuit current, it is evident that the 

2 

decrease of Jg^ due to increasing emitter recombination is about 4.5 mA/cm 
(33.420-28.933=4.49 for N__=10^^ and 30.042-25.547=4.50 for =10^^) . 

liiJ iij3 

The decrease is relatively independent of base lifetime or base recombination 

which is increased by a factor of 100 between these two curves. 

Next, consider the open circuit voltage, V . Column 13 of in 

UU uL# 

Table 2.18 shows that if base recombination dominates over emitter recombina- 
tion, then an increase of emitter recombination from FRAC=10 ^ to 0.8 only 
decreases VOC by 516.85-498. 76=18.1mV, However, if base recombination is 
small or negligible (column 11 of V in the table) , then V is significantly 
reduced by increasing emitter recombination. In this example, the reduction 
is 659. 58-526. 30=133mV or about 5kT/q. 

Two additional effects are illustrated in this calculation, one is the 
ultimate efficiency imposed by interband Auger recombination in the emitter and 
the other is the effect of an acceptor recombination center in the p-type 
emitter instead of the donor just studied. Their effects on efficiency are 
demonstrated in Figure 2.2. 
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Interband Auger recombination in the emitter will slightly lower the 
efficiency when the emitter impurity recombination center density is negligible. 
This is "ilTustrated by curves 1 and lA in Figure 2.2. Curve 1 has no interband 
Auger recombination while curve lA has. Both curves have a very low base 
recombination rate (10 centers/cm ) in order to emphasize the effect of 
interband Auger recombination in the emitter. It is evident that the difference 
between curve 1 and lA at a FRAC (emitter) =10 ^ is only about 1%. This would be 
much less than 1% if higher and' more realistic base recombination is assumed, 
such as 10^^ of curve 11, In the limit of zero emitter impurity recombination, 
emitter Auger recombination would reduce the efficiency by 25.548-23.484=2.06%. 
These are maximum values of curves 1 and lA as FRAC 0 in Figure 2.2. 

The difference between a compensating (donor-like) and noncompensating 
(acceptor-like) recombination center in the boron'dif fused emitter on the 
efficiency is illustrated by the two curves labeled 1-acceptor and 1-donor 
in Figure 2.2. There is negligible difference at low concentration or 'low 
FRAC of these centers. However', when FRAC approaches 1, a donor recombination 
center will compensate the charge of the -diffused boron acceptor impurity in 
the- emitter, thereby reducing the hole density in the emitter and the P+/N 
barrier height.- This will ‘reduce the open circuit voltage which is directly 
related- to the equilibrium potential barrier height of the P+/N junction but • 
will not affect the short circuit current. The solar cell characteristics due 
to -the charge difference of these two recombination- centers in the emitter are 
given in Table 2.19 for the case of FRAC=0.8. The numerical result furhter 
illustrates that ‘VOC is greatly reduced by compensation over that of - non- 
compensation while JSC does not change. 
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Tabel 2.19 Effect of the Charge States of the Recombination Center 
in the Emitter on Silicon P+/N/N+ Solar Cell Performance 
at One AMI Sun. FRAC=0.‘8. 


Computer 

Recombin . 

EFF 

VOC 

JSC 

FF 

Run No. 

Center 

% 

mV 

mA/cm^ 


RUN613 

Acceptor 

15.659 

589.62 

29.101 

.81338 

RUN622 

Donor 

13.752 

527.62 

29.083 

.79689 
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2.6.2 Back Surface Field Layer 

The effect of recombination at the phosphorus -vacancy (PV) center in 
the back surface field layer is to reduce the open circuit voltage from a high 
value determined by base recombination 6r by surface recombination at the 
contact to the value of an equivalent P+/N/<» or P+/N cell with no BSF and with 
infinite 'surface recombination rate. This was already demonstrated in section 
2.2.1 using the gold recombination model. ThePV recombination center does not 
behave too differently. 

A family of curves is again computed to illustrate the reduction of 
efficiency by recombination in the BSF layer and to show how this reduction 
is affected by base recombination. This is given in Figure 2.3. The AMI 
efficiency at one sun is plotted as a function of the fraction (FRAC) of 
phosphorus which is converted into PV pair in the BSF layer. This FRAC is 
also translated into the minority carrier (hole) steady-state lifetime at a 
position 299.20ym which is O.Sym from the back surface and 0.2ym from the 
N/N+ junction. The_ N/N+ junction, is, located--at -IrO- vm -from 'thLe' back surface.' 
Corresponding hole diffusion length at this location is also obtained. 

Both the hole lifetime and diffusion length are indicated' on the FRAC or 
x-axis in Figure 2.3. 

We first note that inclusion of interband Auger recombination in the 
BSF will only reduce the efficiency slightly as indicated by curves labeled 
1 and lA. lA has Auger recombination. This is similar to earlier results. 
However, the compensation effect just described for the emitter layer does not 
appear in the BSF layer to any significant amount. This is due to the fact 
that as FRAC approaches 1, the large density of PV recombination center will 
already have reduced the N/N+ junction to a N/« junction so that compensation 
cannot reduce the effective back surface recombination velocity further. On 
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Figure 2.3 AMI efficiency as a function of phosphorus-vacancy recombination 
center concentration, back surface field lifetime and 
diffusion length in P+/N/N+ silicon solar cells with the 
base recombination center density as the constant parameter. 
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the other hand, if the recombination center in the BSF layer is noncompensating, 
then the N/N+ barrier height would not be reduced as FRAC approaches 1. But 
the recombination rate in the BSF layer would increase to such a high 
value that the N/N+ barrier becomes essentially a N/“> contact. 

Another feature in Figure 2.3 is that the BSF recombination has almost 
no effect on the efficiency when the base lifetime of holes is below about 
lOys which corresponds to a hole diffusion length of about lOOpm. This is 
demonstrated in curve 11. For a hole lifetime in the base of 350ys, BSF 
recombination can reduce the efficiency at most by 1.5% as indicated by 
curve 10. 

The large variation of the steady-state hole lifetime in the base with 

BSF recombination center density, for example from 350ys to 150ys in curve 10, 

comes from the intermediate injection level condition in the base. At low 

FRAC or low base recombination, the injection level is high so that the hole 

lifetime approaches the high level lifetime which is T^Q+Tp^=361+50=411v]s 

10 3 

from Table 2.3 for N„„„=10 centers/cm . As FRAC increases toward one, 

lib 

the BSF recombination also increases so that the injection level in the base 

decreases due to the lower hole concentration reaching the N/N+ high-low 

junction. This would decrease the hole lifetime toward the low injection 

level value which is t ,=50ys. 

pi 

The reduction of the efficiency with increasing BSF recombination or 
FRAC, indicated in Figure 2.3, arises mainly from a decrease of VOC just 

stated. This is further demonstrated by the numerical results tabulated in 

10 3 

Table 2.20 for N„„„=10 centers/cm , 

lib 
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- Table 2.20 Effect of the PV Acceptor Recombination in the Back Surface 
Field Layer on Silicon P+/N/N+ Solar Cell Performance at 
One AMI Sun. 


Computer 
Run No . 

FRAC 

EFF 

% 

VOC 

mV 

JSC 

mA/ cm^ 

FF 

RUN671 

10"^ 

16.958 

638.27 

32.650 

.72359 

RUN672 

10"'^ 

16.946 

632.01 

32.646 

.73035 


10“^ 

16.860 

611.98 

32.619 

.75099 

RUN674 

10"^ 

16.542 

583.27 

32.461 

.77810 

RUN675 

10”^ 

15.973 

560.12 

32.037 

.79151 

RUN680 

0.99 

15.664 

551.71 

31.748 

.79522 


The tabulated results show that recombination in the BSF layer significantly 

reduces VOC, from 638.27 to 560.12 or 78.2mV, but the reduction of JSC is 

2 

rather small, from 32.650 to 32.037 or 0.61mA/ cm when the PV recombination 
center density is increased to 10% of that of the phosphorus concentration in 
the diffused BSF layer. The last row for RUN680 demonstrates the absence of the 
performance dropoff when 99% of the phosphorus donors in the BSF layer are 
converted into phosphorus-vacancy acceptor recombination centers . This is 
in contrast to the large compensation effect by the boron-vacancy donor 
complex on the boron acceptor in the diffused emitter discussed in the 
preceding subsection and shown in Figure 2.2. 

A useful guide can be deduced from the results in this table. The 
improvement of efficiency and VOC that can be achieved by reducing the PV 

2 

recombination center in the BSF layer is not very large. Based on 95 mW/cm 

2 

AMI incident power, efficiencies of 16.8%, with VOC=610mV and JSC=32. 2mA/cm , 
have been obtained from Sandia's high efficiency cells. ^ Comparing these nianbers 
with those in Table 2.20 suggests that the concentration of the electrically 
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active PV centers which, are also active recombination sites in the phosphorus 

_3 

diffused BSF layer cannot be higher than about 10 times the phosphorus 

concentration. This is substantially smaller than the value, PRA.C=0.6, that 

2 

is assumed in the detailed analysis of the phosphorus diffusion profile. 

In order to bring these into agreement, the capture rates of electrons and 
holes at the PV center given in Table 2.3 must be reduced by 600. This would 
reduce the capture radius of holes by the negatively charged PV center from 
25A assumed earlier to lA. The value of lA is inconsistent with the physical 
size of the PV center which must occupy two silicon lattice sites or about 5A. 
It is also inconsistent with the long range Coulomb force of the negatively 
charged PV pair. Alternatively, the active PV recombination center density 
may be indeed very low due to slow cooling and other annealing treatment 
after the lOOOC phosphorus diffusion. 
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3. 


TECHNICAL DISCUSSION - Ti -DOPED CELLS 


3.0 Objectives 

Substantial efforts have been spent to study the effects of metallic 
impurities on the performance of silicon solar cells since they are not 
completely removed in solar cell grade silicon single crystals due to cost 
considerations. These studies were carried out by Westinghouse and previously 
by Monsanto groups. Recently, accurate minority carrier diffusion length 
and lifetime measurements in the base of solar cells made on impurity-doped 
silicon have been made by the Northrop group. The crucial parameters needed 
to characterize the recombination properties at these metallic impurities are 
the thermal capture and emission rates of electrons and holes. Alternatively, 
it is sufficient to completely characterize the recombination properties at 
an impurity recombination center by having three parameters : the thermal 
capture rates of electrons and holes and the thermal activation energy. But, 
these parameters are still not accurately determined. Measurements by the 
so-called DLTS (Deep Level Transient Spectroscopy) method have given qualitative 
indications of only the energy levels. Accurate values of these parameters can 
only be obtained by the t/uin&^znt capacXtonaz m&thocU which employ the 
following techniques ; 

TSCAP - Thermal Stimulated Capacitance 

VSCAP - Voltage Stimulated Capacitance 

LSCAP - Light Stimulated Capacitance 

These methods will provide all of the recombination parameters, .including both 
the thermal emission rates, activation energies, and -the. t-ho/rndt dCLptuAQ. Kojtu 
aJ: thoMmaJi dqiuJLLbnMm. 

Without these parameters, it is difficult to accurately predict the 
maximum allowable concentration of a metallic recombination impurity to give a 
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prescribed solar cell performance. However, a nmber of experimental measure- 
ments have been made on Ti-doped silicon solar cells so that an extrapolation 
can be made to give a preliminary set of values for these recombination 
parameters at the Ti recombination centers. These measurements include the 
Ti concentration in each crystal, the minority carrier diffusion length in 
the base from which the minority carrier lifetime can be computed, the 
extrapolated minority carrier base lifetime using an experimental curve of 
open circuit voltage decay lifetime versus diffusion length lifetime 

(t^) , the concentration of carbon and oxygen in the starting silicon crystal, 
and some preliminary but fairly reliable thermal emission rates and energies 
obtained by us in Ti-doped p-type Westinghouse silicon single crystals . 

In this chapter, a study of the Ti on the silicon solar cell performance 
is undertaken using these preliminary recombination parameters. The study is 
aimed to provide an indication of the range of solar cell performance one 
can expect in the presence of Ti impurity. It also serves to indicate the 
type of projection that can be made from the computer model and the ambiguities 
which need to be resolved by the transient capacitance measurements . 
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3. 1 Approach 

The approach undertaken here is similar to that described in section 

2.1 of the previous chapter for high efficiency cells. The same computer 
model will be used but the Ti recombination model will be used instead 

of the gold recombination model used previously. The details of the Ti 
recombination model will be described in subsection 3.1.1. The device and 
material parameters to be used are those- from the processing data sheet of 
Westinghouse so that the computed solar cell performance characteristics can 
be compared with the experimental data given by Westinghouse. Since the 
experimental cells have no anti— reflection coating, the surface reflection 
loss at the silicon surface is included in the photogenerators of the model. 

3.1.1 Ti Recombination Model 

One electron and two hole thermal capture coefficients at the Ti 

centers in Ti-doped silicon are obtained from the base lifetime versus Ti 

concentration data. These data are shown in Figure 3.1. The Ti - concentration 

21 

data are taken from the Westinghouse Quarterly Reports and the Monsanto Pinal 

22 23 

Report. The base lifetime data are taken from the Northrop Final Report 

which were computed by them from their diffusion length measurements on cells 

made on the Tl-doped silicon of Monsanto and Westinghouse, For the more 

recent Westinghouse silicon doped with Ti which had no diffusion length 

measurements made by Northrop (crystal number greater than W056), the curve 

of (Westinghouse data) versus r (Northrop diffusion length data) given by 

Northrop is used to compute the base lifetime since for each cell is given 

by Westinghouse so that the base lifetime, can be read off the vs 

curve. This procedure provides an estimated base lifetime for all of 'the more 

recent Westinghouse cells. 

The points in Figure 3.1 can be separated into three groups. The lowest 
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Steady-State Base Lifetime (sec) 









line contains points' from p-base silicon cell's doped with Ti. The ' lifetimes 
correspond to electrons (minority carriers) captured into the E^+336 mV level. 
This energy level and the other two at E -470 mV 'and E -182 mV were determined 

■ L* 

by VSCAP transient capacitance measurements on Ti-doped silicon using Mg/p-Si 

Schottky barrier diodes. The assignment of these data to this 

transition is consistent with the position of the bulk Fermi level in the p-type 

silicon which is E^-E^ = 168-^229 mV at 297K. Thus, the Fermi level position 

is below this Ti level which is located at = 336 mV. Since the lifetime 

is small, it is also consistent to assume that these are low level lifetimes in 

the base. In addition, it is not inconsistent with some of the preliminary 

experimental indications of Westinghouse that emitter recombination is probably 

not too important even if the N+ phosphorus diffused emitter has some gettering 

action of Ti but this point must be more carefully investigated. 

The middle and the highest lines on Figure 3.1 contain two points each 

from Ti-doped n-type silicon, one from Westinghouse and the other from Monsanto. 

Although there are two points each, they are unmistakenly separated even allowing 

the large uncertainty of the data which we have not shown in the figure. In 

addition, there are many points from p-type Ti-doped Westinghouse crystals 

which are clustered around the middle line. These are discussed later. 

For the four points from n-type silicon, a distinction among them is 

that the two crystals on the middle line contain either a high concentration of 

carbon or a low concentration of oxygen, while the two crystals on the upper 

line have just the opposite proportional of carbon and oxygen. Their carbon 

21 22 

and oxygen concentrations were given by Westinghouse and Monsanto and are 
given in Table 3.1 on the next page. 

The basis for these two groupings to give the two upper lines in 
Figure 3.1 is as follows. The upper line for the two crystals, W065 and C17B, 
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Table 3.1 Carbon and Oxygen Concentrafc-ion in-.n-Si 
Doped with Ti. 


Crystal 

Number 

Resistivity 

ohm-cms 

Dopant 
lOf^P/cm? • 




W0&5N 

1.8 

3. 

2 or 3.5 

<2 

176 

W121N 

1.8 

3. 

0.4 

4.7 

121 ’ 

C17B 

0.6 

10. - 

9.5 

1.5 

174 ’ - 

. M17B 

0.54 

11. - • 

<1.3 ■ 

■0.-6. 

5.6 - 


were grown by the Czochralski method and have high oxygen concentration (about 

1.8x10 O/cm ). They also have low carbon concentration. On the other hand, 

for the two crystals on the middle line, crystal W121 has a relatively higher 

» - > 

carbon concentration and somewhat lower oxygen concentration, while crystal 

M17B has low concentration of both oxygen and carbon. The following tentative 

model seems to emerge. The two high oxygen crystals with low carbon content, 

W65 and C17B, are likely to have a shallow acceptor level associated with the 

oxygen which is present. Such a shallow level, in the form of an oxygen- 

vacancy pair, is known to form when vacancies are generated in oxygen-rich 

silicon crystals. Vacancy is very mobile and gets trapped next to the oxygen. 

The high Ti concentration during crystal growth may creat a high density of 

vacancies which may be trapped at the oxygen sites upon cooling. This model 

is suggested by the well-known oxygen- vacancy energy level in silicon which 

is E -174 mV and it is nearly equal to observed activation energy of an acceptor 
u 

level, E — 182 mV, using the VSCAP method. Thus, this upper level and line 

can be attributed to the OV center but a more complex 0 V Ti center cannot be 

X y z 

ruled out . 

The association of the E -182 mV level with the OV pair or oxygen is 
also consistent with the two crystals lying on the middle line, W121 and M17B. 
The M17B crystal, grown by float-zone method, has a low concentration of both 
oxygen and carbon. Thus, the oxygen-vacancy acceptor will have very low 
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concentration, if it exists in M17B. The hole lifetime in the n-ty.pe base cell 
of this crystal is then associated with the hole capture into the middle 
level, E„“470 mV since the E -182 mV level has low concentration. 

Ci 0 

For the W121 crystal, the carbon concentration is high, although not 

substantially higher than the other three crystals. But the higher carbon 

concentration suggests that the OV level (E -182 mV) is impeded by the 

o 

presence of carbon. For example, -carbon may., form a complex with the OV pair 
and either render it electrically inactive or drastically change its energy 
level. Thus, the hole capture transition is again dominated by the E^-470 mV 
level in this crystal due to the lack of the E^-182 mV level. 

There are four Ti-doped p-Si from Westlnghouse, W102, W122, W123 and 
W140, whose points are clustered around the middle line for the E^-470 mV level. 
At least two possibilities exist to account for these points; (i) the E.^+336 mV 
level does not exist in these crystals and the transition involves electron 
capture into the E -470 mV level, and (ii) high level injection level occurs 
in the base of the cells from these crystals so that hole capture into the 


mV level dominates the recombination. 


It is obvious that all of these considerations can be readily checked if 


transient capacitance measurements are made on these crystals. 

The capture rates can be readily computed from the slope of the straight 
lines of Figure 3.1. The corresponding thermal emission rates for the majority 
carriers were obtained from VSCAP measurements. These results are listed in 


Table 3.2 on the next page and are used for further calculation to give 
the complete set of recombination parameters at thermal equilibrium. 

The majority carrier capture rates can be estimated by assuming a 
temperature independent capture cross section and by assuming that the emission 
rates are Independent of the electric field. Then, e^=c^N^exp(-E/kT) and 


- 67 - 



Table 3.2 Experimental Thermal Emission and Capture Rates 
in Ti-Doped Silicon 


E^ + 336 

mV 

Level 


c 


3.7x10"^ 

3/ 

cm /sec 

n 




e 

P 

= 

1.1x10^^ (T/300) ^exp (-336mV/kT) 

1 /sec 

\ - "0 

mV 

Level 


e 

n 

- 

5 . 6x10^^ (T/300) ^exp (-470mV/kT) 

1 /sec 

c 


2.5x10"® ■■ 

3/ 

cm /sec 

P 




\ 

mV 

Level 


e 

n 

= 

3 . 6x10^^ (T/ 300) ^exp (-182mV/kT) 

1 '/sec 

c 

— 

-9 

5. xlO 

3, 

cm / sec 


e =c N„exp(-E/kT) . These can be nsed with the experimental emission data to 
P P V 

compute the thermal capture rates, c^ and c^. The fourth parameter, the 

minority carrier thermal capture rate, of each of the three levels is computed 

2 

using the mass action law, e e =c c n., where three of the rates are already 

n p n p X •’ 

obtained and n^ is also given. 

The solar cell calculations are made at 297K, The kinetic recombination 
coefficients computed using the procedure just described are obtained for the 
three levels observed in Ti— doped silicon solar cells. These are given in 
Table 3.3 on the next page. 

The choice of the charge state shown in Table 3.3 is based on the 
relative size of the electron and hole capture rates and the position of the 
level to a band edge. The two lower levels are assumed to be donor like since 
c^>>Cp for both levels. The upper level is assumed to be acceptor like since 
it is closer to the conduction band edge and since the OV center is an acceptor. 
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Table 3.3 Derived Recombination Parameters of the Ti Related 
Energy Levels in Si at 297K. n^=l. 036x10 ^ ^/cia^". 


E^-182 Acceptor, Level 


14 3 ' 

Lifetime at 10. Levels /cm 


c 

— 

,13 

E-07 

3 / 

cm /sec 

T 


760 

ns 

n 




n 



e 



E+09 

1 /sec 





n 








c- 


.5 

E-08 

3-, 

cm’ /sec 

-T 


2000 

ns- 

P 




P 




e 


..22 

E-04 

1 /sec 





P 








E -470 Donor 

Level 







t; 







- - 


c 



.,2 

E-05 

3, 

cm / sec 

T 

_ 

5 

ns 

n 




n 




e 

= 

. 6 

E+06. 

1 /sec 





n 








c 

- 

..25 

E-07 

cm /sec 

T 

= 

.400 

ns 

P 




P 




e 

= 

.92 

E+01 

1/sec 





P 








E^+336 Donor 

Level 







c 


.4 

E-06 

3 / 

cm /sec 

T ■ 

— 

25 

ns 

n 




n 




e 

= 

.16 

E+01 

1 / sec 





n 








c 


..82 

E-08 

3 , 

cm /sec 

T 


1150 

ns 

P 




P 




e 

= 

.22 

E+06 

1 /sec 






P 


In the computer calculation, we shall use the two lower levels and 
assume that they are coupled and come from one species of Ti recombination 
centers. No experimental data are available to show if the two lower levels 
are associated with a single Ti recombination center or they are two independent 
energy levels associated with two different species of Ti recombination centers. 
For example, in the two independent level model, one of the two Ti species can 
be an isolated substitutional Ti center while the other, a Ti-0 or Ti-A pair 
where A is an unknown impurity. 
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3.1.2 Computer Model 

The computer .mpdel is nearly identical to that used for the high 
efficiency cell studied in chapter 2 and described in section 2.1.2. The only 
difference is that the length of the transmission line is 250 pm instead of 
300 pm • since the experimental Westinghouse cells have a nominal thickness of 
250 pm. For the p-base cell, the emitter is 0.25 pm and for the n-base cell, 
the emitter is 0.185 pm. There is no BSF layer in either of these cells. 
However, the four region model is still used in which the fourth region, 
occupied by the BSF layer previously, is just an extension of the third, 
region or the quasi-neutral base region. In projecting the maximum allowable 
Ti concentration for a desired efficiency, cells with BSF are used in the 
computer model since they give the best efficiency with the highest open 
circuit voltages. For these BSF cells, the four region model is needed and 
the BSF layer is assumed to be 1 pm thick while the fourth region is 2 pm thick. 

The better integration scheme, exact for dx and Simpson for dX, is used 
for all of the calculations made for Ti-doped cells. 

•3.1.3 Device and Material Parameters 

The parameters in Table 3.4 are used so that results of this study 
can be compared with the experimental data of the Westinghouse Ti-doped cells. 
The parameter values are obtained from the Quarterly Reports of Westinghouse. 

The fabrication history as well as the electrical properties of Che finished 
silicon solar cells, doped with metal impurities are extremely well documented 
in these quarter reports. 

Three columns are given in Table 3.4, one for the p-base cell, one for 
the n-base cell and the third for the highest efficiency Ti-doped P+/N/N-f- cell. 
The separate listings of the experimental cells is due to the slight differences 
in the base dopant impurity concentration, surface concentration and profiles.' 
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Table 3.4 Device and Material Parameters of Ti~Dop,ed 
Silicon Solar Cells of Westingbouse. 


Item 

. Parameter 

— 

p-Base 

n-Base 

P+/N/N'f 

1 

T 

(K) 

297 

297 

297 

2 

n . 

X 

(cm~^) 

0.1036E+11 

0.1036E+11 

0.1036E+il 

3 

. L 

(ym) 

250 

■ 250 

250 ■ 

4 


(ym) 

■0.250= 

0.185 

0.250 

5 

^BSF 

(ym) 

- 

- 

- 

6 


(ohm-cra) 

. 3.5 4.5 

1.4 2.2 


7 


(cm~^) 

5.0E-f20 

■ 2.5E+20 

5.0E-+-20 

8 


(cm 

5.0E+15 

5 . OE-115 

-5.0E-fl5 

9 


(cm”^) 

- 

- 

2..^E+20 

10 

Emitter 

Profile 

1-X (two) 

1 - 


11 

BSF Profil'e 

... 

- 

- 

exp(-*X^) 

1 

'l2 

Npi 

(Ti/cm^) 

constant 

constant 





4E11 ^ 4E14 

4.0E+13 



For the phosphorus diffused emitter of the p-base cell, the diffusion . 

g 

profile was obtained by sheet resistance- anodic etching method. It contains two 

linear segments due to the particular diffusion schedule and ambient gas used. 

This profile drops off slower than- the exp(-X^) profile used previously for the 

phosphorus diffused N+ layer. The boron diffused emitter profile was also obtained 

9 

by sheet resistance-anodic etching method. It is essentially identical to the 
2/3 

1-X profile used earlier for the high efficiency cell studies in chapter 2. 
Slight differences between these empirically fitted profiles and the experimental 
profiles will not significantly affect the final results and conclusions as we 
have demonstrated in the high efficiency cell studies in Chapter 2 if the 
enhancement of Ti solubility in the emitter is not large which is apparently the 
case for low Ti concentrations. For high Ti concentrations (4x10 Ti/cm ), 
enhancement of solubility in the emitter may be large and the solar cell 
characteristics will be' more critically dependent on the emitter diffusion profile. 
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3.1.4 Solar and Optical Properties 


The AMI solar spectra, silicon absorption coefficient, C6(l), and silicon 
surface reflection coefficient, R(A), used in this study are those given in 
Table 2.7 of section 2.1.4. To- simulate the performance of the Westinghouse 
cells,' silicon surface reflection loss is taken into account by reducing the 
incident solar spectral irfadiance [the column PC^) in Table 2.7] by the factor 
l-R(X). This is done by including the factor l-R(X) in the integrand of the 
distributed photocurrent generator given by Equation (2.1). 

•The maximum AMI phbtocurrent that can be collected in a 250 ym thick 

silicon solar cell with no anti*-ref lection coating is calculated to be 22.025 

2 2 
mA/cm . Without surface reflection loss, it is 34.198 mA/cm . Thus, the 

surface reflection loss at a bare silicon surface amounts to a reduction of 

[1 - (22.025/34.198.)] = 0.356 or 36%. However, this reduction factor is not 

a constant since it will depend on the cell thickness as well as the distribution 

of the recombination centers within the cell because of the wave length 

•dependence of the absorption and reflection coefficients. But, the figure of 

0.36 or 36% can be used as a. rule of thumb estimate when considering the 

short circuit current and the efficiency. An estimate of the increase of the 
< . ' • * 

open circuit voltage due to the higher photocurrent when there is no surface 
reflection loss can also be made. If we assume low level injection condition 
inside the cell, then 

AV_ = (kT/q)log (34.198/22.025) = 11.37 mV 
For high injection level, the increase from an uncoated to a perfectly coated 
cell would be twice this value or 22.75 mV. 
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3.2 Results on p-Base Cells 

In order to compare the computed results with the experimental -data of 
the Ti-doped silicon cells fabricated by Westinghouse, the- diffusion profile 
and the dopant concentration are selected to match their reported data. From 
Table 3.4, the net ionized impurity concentration profile used in the computer 
program is given by 

^DD “ ^AA " 5x10^'^ [1.04 - 3.29(Z/L^).] - 5x10^^ ' 0 < Z <■ 0.07694pm- ' 

■= 2x10^^ [1.00 - 1.00(Z/L^)] - 5x10^^ 0.07694 < ,Z < 

- -5x10^^ ■ ■ L^< Z < 250pm (3.1) 

Here, [l-25Cg/CQ]=0’. 2500626156 pm and the large number of significant 

figures given is to indicate that there is no jump in the N -N value at 

L)JJ AA 

the junction Xj, only a change in slope or concentration gradient. This 
impurity concentratiion profile fits the experimental profile obtained. by • 
Westinghouse for the phosphorus diffused emitter very well as demonstrated- in 
Figure 3.2. • ' . 

The Ti concentration is assumed spatially constant but as we shall see 
from the results , there is an indication of Ti pileup in the phosphorus 
diffused emitter for the highest Ti concentration cell from Westinghouse 
crystal WOOS. 

The cells are not coated with anti-reflection coating, thus, the surface 
reflection loss is included in the calculation by the procedure described in- 
section 3.1.4. 

The computed results for the N+/P cells are compared with Westinghouse 's 
21 ■ ■ 

experimental data. Comparisons are made for EFF, VOC, JSC, and FF in 

Figure 3.3(a) to (d) as a function of Ti concentration. Data of only four 

/ 

Ti-doped N+/F cells are included which are from Westinghouse crystals 
W048, W033, W066 and W008, in increasing Ti concentration. 
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Phosphorus Concentration P/cm 
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0.4 
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1.0 


Relative Position from Surface ’X/X 

J 


Figure 3.2 The phosphorus diffusion profile (experimental points) 
and the two constant concentration— gradient empirical 
fits (solid curve) to experiments. X =0,25ymt 
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The best experimental- theoretical, agreement is obtained for the open 
circuit voltage, shown in Figure 3.3(b). This is probably due to its insensi- 
tivity to random experimental variations of device and material parameters 
since it is measured at zero current. The agreement also indicates self- 
consistency in the assumption that base recombination dominates in these 
Ti-doped cells. 

The poorest agreement is in the fill-factor. This is. probably due to 
its high sensitivity to the contact resistance at both the front and back 
surface contacts and also its dependences on the high conductance paths which 
may have a high occurrence rate in the emitter of cells highly doped with Ti. 

The efficiency and the short circuit current show better theoretical- 

experimental agreement at low Ti concentrations as revealed in Figures 3.3(a) 

14 

and (c) . For the highest Ti concentration cells (crystal WOOS, N^^=3.8xl0 
3 

Ti/cm ), the experimental values are more than 50% below the theory, although 
VOC [Figures 3.3(b)] and FF [Figure 3.3(d)] are nearly equal to the theory.' 
Thus 5 the smaller experimental EFF and JSC can only be due to some additional 
recombination mechanisms not included in the computer model. This contention 
is further supported by the nearly ideal dark I-V slope of one, I exp (qV/kT) , 
over several decades of current, including the open— circuit voltage and short- 
circuit current points. One possible cause is the high emitter recombination 
current due to enhanced solubility of Ti in the high concentration phosphorus 
diffused emitter. This would be more important when the silicon is doped to a 
high Ti concentration since then a source of high Ti concentration is present. 
Enhancement of Ti solubility in the emitter was not included in the computer 
model used to compute the theoretical curves given in Figure 3.3. However, 
sample calculations were made which showed that even with an emitter lifetime 
of T (0.2iim)=4 ps at a distance of 0.05pm from the junction, JSC is reduced 
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2 « ' " 2 

only ^to 10 mA/cra from' 12.5 mA/'cm , the latter given in Figure 3.3(c) ‘at a 

14 3 

bulk Ti concentration of 3.8x10 Ti/cm . This is still 25% higher than the 
= 2 

experimental value of 7.7 mA/cm . Thus, increased emitter recombination' 
current due to enhanced Ti solubility in the emitter appears to account for 
only half of the theoretical-experimental difference. Another possible 
source is areal inhomogeneity of the Ti concentration. This could make 
part of the cell area substantially higher in Ti -concentration and lower in 
short circuit current and efficiency so that the JOC and EFF of the entire cell 
are considerably below that predicted using an average Ti concentration. 
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Figure 3.3 Comparison of calculated and experimental characteristics 
of Ti-doped p-base silicon N+/P solar cells at one AMI sun. 

(a) Efficiency, EFF 

(b) Open-circuit voltage, VOC 
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Figure 3.3 Comparison of calculated and experimental characteristics 
of Ti-doped p-base silicon N+/P solar cells at one AMI sun. 

(c) Short-circuit current, JSC 

(d) Fill-factor, FF 


3.3 


Results on n-Base Cells 


Very few n-type silicon doped with Ti have been grown thus far. 
Extensive data on cell performance have been documented only on two crystals, 
W065N and W121N. A comparison of theory with experiment will be made only 
with cells made from crystal W121N where the middle Ti donor center is 
believed to dominate the recombination current in the base. Computations 
have not been made for cells made on W065N since it was thought that the 
major base recombination current in this crystal is via the oxygen-vacancy 
or oxygen related recombination center. 

The net ionized impurity concentration profile, using Table 3.4, is 
given by 

Naa - Npp = 2.5xl0^°[l - (Z/L^)^^^] - 5x10^^ (3.2) 

The Ti concentration is again assumed to be spatially constant and the silicon 
surface reflection loss is included in the calculation. The theoretical- 
experimental comparison is given in the following table. 

Table 3.5 Comparison of Experimental and Computed Performance of a 

Ti-Doped P+/N Si Solar Cell at One AMI Sun and 4x10^^ Ti/cm^ . 


1 

Cell No. 

EFF 

% 

VOC 

mV 

JSC 

mA/cm^ 

FF 

W121N 

Expmt 

7.51 

525 

19.9 

.681 

RUN733 

Theory 

8.54 

536 

17.6 

.804 


The agreement is fairly good considering that several experimental 
parameters are not accurately known and the light source used is not identical 
to the AMI solar spectra. The largest difference is in the fill factor, FF. 
The low experimental value is probably again due to contact resistances. The 
slightly higher experimental JSC probably reflects the difference between the 
light source and the solar spectra. 
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3.4 Projected Maximum Performance 

In view of the fairly good agreement between experimental data and 
theoretical prediction presented in section 3.2 and 3.3 for both the p-base 
and the n-base Ti-doped cells, the Ti recombination parameters given in 
Table 3.3 are used to compute a family of theoretical performance curves as 
a function of the Ti concentration. Since the hole capture rate (2.5E-8) 
at the upper Ti donor level (E^-470mV) is smaller than the electron capture 
rate (4E-7) at the lower Ti donor level (E^+336mV) , the efficiency is expected 
to be higher for the n-base than the p-base cell. The highest efficiency 
cell would then be the P+/N/tW- geometry with back surface field and anti- 
reflection coating. Thus, two theoretical curves are computed, one representing 
the worse case of a N+/P cell without anti-reflection coating and the other, 
the best case of a AR/P+/N/N+ BSF cell with perfect anti-reflection coating. 

The efficiency of these two cases as a function of Ti concentration is given 
in Figure 3.4 at the one AMI sun condition. The cell geometry and diffusion 
profile used are tabulated in the third column of Table 3.4. 

The results of Figure 3.4 show that for an efficiency of 16%, the 
maximum allowed Ti concentration is about 6x10 Ti/cm in the most efficient 

cell geometry, i.e. the P+/N/N+ cell with anti-reflection coating. The other 

2 

characteristics of this cell are: Vq^= 576 mV, Jg^=30.2 mA/cm and FF=0.8215. 

A single point is also computed for a P+/N cell without anti-reflection 

coating and shown in Figure 3.4. Table 3.6 gives all the cell characteristics 

13 3 

for the three cells at 4x10 Ti/cm concentration. The differences between 
the more efficient P+/N cell and the less efficient N+/P cell come from 
recombination in the base region. In the n-base cell of RUN731, the hole 
capture rate at the upper donor level, E^-470mV, is smaller while in the 
p-base cell of RUN732, the electron capture rate at the lower donor level. 
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Figure 3.4 Theoretical AMI efficiency as a function of Ti 

concentration for bare surface N+/P, bare surface 
P+/N and anti-reflection coated P+/N/N+ silicon 
solar cells at 297K. 
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Table 3.6 Performance Comparison of Bare N+/P, Bare P+/N and AR 
Coated P+/N/N+ Silicon Solar cells At One AMI Sun and 
Containing 4xlO^^Ti/cm^. 


Computer 


JPHOTO 

EFF 

VOC 

JSC 

FF 

Run No. 

Cell Type 

mA/ cm^ 

% 

mV 

mA/ cm^ 


RUN729 

AR/P+/N/N+ 

34.198 

13.691 

547.85 

27.479 

.80869 

RUN731 

P+/N 

22.025 

8.552 

536.12 

17.630 

.80462 

RUN732 

N+/P 

22.025 

6.708 

482.49 

15.595 

.79271 


is larger. Thus, the base lifetime in the n-base cell is larger than that 


in the p-base cell. High injection level effect will make the difference 

“7 —8 

in lifetime smaller than the capture rate ratio of c /c =3.7x10 /2.5xl0 

n p 

=14.8, using the values from Table 3.3. 
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4. 


CONCLUSION 


Interband Auger recombination in the diffused emitter and back surface 

field layers does not affect the silicon solar cell j)erformance significantly 

even if the Auger recombination rate data are in error and too small by a factor 

of 100. It becomes important only when the carrier concentration in these 

20 3 

layers is higher than about 10 carriers/cm throughout the entire layer. 

Surface recombination at the front and back surfaces has negligible 

effect .on cell performance unless the interface recombination velocity at 

these surfaces lies between the diffusion velocity, 10^ cm/s, and the effective 

3 

recombination velocity, 10 cm/s, of the emitter and back surface field layers. 

Enhanced solubility of the metallic recombination impurities in the 
emitter and the back surface field layers due to high phosphorus and boron 
concentrations in these layers can reduce the short circuit current and the 
open circuit voltage respectively if the solubility reaches 10% of the diffusion 
impurity (boron and phosphorus) concentrations. 

There is little difference in solar cell performance computed from using 
the Gaussian diffusion profile and the measured phosphorus and boron diffusion 
profiles. 

For the Sandia high efficiency (16.8% AMI) and high resistivity (10 J^-cm) 
P+/N/N+ silicon cells, 1% of gold in the diffused emitter or 0.1% of phosphorus- 
vacancy or E-center in the diffused back surface field layer will account for the 
observed efficiency, open circuit voltage, short circuit current and fill factor. 
The theoretical maximum AMI efficiency for this cell design is 19.26% when 
recombination in the emitter and back surface field layers is negligible. 

However, if recombination in the BSF layer is the limiting factor, then eliminating 
it entirely will only improve the efficiency by 0.1% to 16.9% and the open circuit 
voltage by 26 mV to 638 mV because high injection level condition would prevail 
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in the entire base which would reduce the fill factor. On the other hand, if 

emitter recombination is the limiting factor, then eliminating it entirely 

would improve the efficiency to 19% and open circuit voltage to 660 mV. 

The reported AMI efficiency of 16,8% by Sandla is based on measurements 

using a solar simulator while the computed value is based on the measured solar 

spectra. Thus, some upward correction of the reported efficiency may be needed. 

This brings the reported efficiency to almost the theoretical maximum, 19%, 

of this cell design. Large improvement of performance cannot be- expected in 

this cell design which already has a base lifetime of 700 ys. 

Of the three energy levels detected in Ti-doped silicon and silicon 

solar cells, the* two lower levels are donor like and related to the Ti impurity 

while the upper level is acceptor like and probably related to the presence of 

oxygen and carbon. The maximum. allowable Ti concentration for a 16% AMI 

12 3 

efficiency is about 6x10 Ti/cm in the best cell type, P+/N/N+, which is 

the least sensitive to Ti contamination. A cell thickness of 250 ym and a 

15 3 

base doping of 5x10 pho‘sphorus/cm or a resistivity of 1.5 ohm-cm are also 

assumed. This cell would have an open circuit voltage of 576 mV and a short 

2 - 

circuit current of 30 mA/cm . 
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